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ABSTRACT 


The  specific  nc  flow  resistance,  characteristic  impedance  and  oomplex  wave 
number  of  air-filled  fiberglass  have  been  extracted  from  measurements  of  the 
peak  and  half-power  frequencies  of  the  two  lowest- frequency  normal  modes  of  a 
rectangular  plexiglass  cavity  partially  filled  with  fiberglass.  Using  Biot 
theory  [Ref's  8,  9  and  101  to  describe  the  sound  field  within  the  fiberqlass, 
a  dispersion  relation  for  the  allowed  modes  in  the  cavity  was  derived,  from 
which  the  above  properties  were  extracted.  The  values  determined  by  this  new 
method  compared  favorably  (within  25%)  with  the  results  of  measurements  using 
classical  methods.  The  new  method  differs  from  the  classical  resonance  tube 
method  [Ref.  4]  by  1)  the  use  of  a  wave  theory  to  descibe  the  sound  field 
within  the  fiberqlass  and  2)  the  use  of  a  mode  with  particle  velocity  parallel 
to  the  air-fiberglass  interface.  The  differences  make  it  possible  to  extract 
the  acoustic  properties  of  bulk  fiberglass  using  only  a  thin  sample. 
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I .  INTRODUCTION 


A.  BACKGROUND 

Elastic  wave  propoqation  in  a  fluid-saturated  porous  solid  is  a  very 
important  and  active  subject  of  research  with  a  wide  spectrum  of  applications 
including  oil  prospectinq,  nuclear  test  ban  treaty  verification  and  sound 
transmission  and  reflection  at  the  ocean  floor. 

When  the  fluid  in  the  pores  is  highly  oompressible  (air  for  example)  the 
acoustic  properties  are  completely  determined  by  the  characteristic  impedance, 
Zchar'  or'  QQuivalently  ty  the  complex  wave  number,  k  .  From  these  quantities 
can  be  extracted  parameters  which  depend  only  on  the  microscopic  geometry  of 
the  porous  solid,  the  values  of  which  can  then  be  used  when  the  fluid  in  the 
pores  is  not  hiqhly  compressible  (water  for  example).  It  is  therefore 
desirable  to  have  accurate  and  reliable  methods  by  which  these  parameters  can 
be  measured. 

B.  HISTORY  OF  CHARACTERISTIC  IMPEDANCE  AND  COMPLEX  WAVE  NUMBER  MEASURING 

TECHNIQUES 

Rayleiqh  (1929)  [Ref.  1,  p.  328]  has  investigated  the  performance  of  an 
idealized  absorber  of  sound  consisting  of  snail  bundles  of  tubes.  He  extended 
the  theories  of  Helmholtz  and  Kirchoff  that  relate  to  the  transmission  of 
sound  alonq  tubes  and  was  able  to  deduce  a  value  for  the  absorption 

coefficient  at  the  surface  of  such  a  material  in  terms  of  the  propagation  of 

the  surface  area  occupied  buy  the  open  pores,  the  radius  of  the  tubes,  and  the 

physical  constants  of  air.  This  analysis  has  provided  a  basis  on  which  to 
found  more  general  theories  and  description  of  sound  propagation  through 

porous  materials. 
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Scott  (1945)  (Ref.  21  proposed  a  method  of  determining  the  complex  wave 
number,  k  which  describes  the  spatial  variation  of  the  acoustic  disturbance  in 
a  porous  medium.  This  method  consists  of  measuring  the  attenuation  constant 
and  the  wavelenqth  constant  of  a  plane  progressive  wave  by  probing  the  porous 
medium  with  a  probe  tube  microphone.  Scott  then  determined  the  chatacteristic 
impedance  of  the  material  with  an  apparatus  devised  by  H.  0.  Taylor  (1913) 
(Ref.  31  to  measure  the  sound  damping  effects  of  acoustical  materials.  This 
apparatus  is  similiar  to  the  present  day  impedance  tube. 

Beranek  (1949)  (Ref.  4]  describes  several  methods  of  obtaining  a 
material's  acoustic  impedance  by  measuring  the  properties  of  the  sound  field 
in  the  vicinity  of  an  absorbing  material  exposed  to  a  normally  incident  sound 
wave.  Two  of  the  most  tractable  methods  are  based  on  transmission  line 
theory:  the  impedance  tube  method  and  the  resonance  method  (normal 
incidence),  with  the  latter  method  generally  being  viable  only  for  determining 
surface  impedance.  Beranek  also  developed  a  method  to  calculate  the  complex 
prooaqation  constant  from  measurements  of  a  material's  characteristic 
properties,  i.e.  porosity,  pore  size,  etc.  Once  these  properties  are 
determined,  the  characteristic  impedance  can  be  calculated. 

More  recently,  Yaniv  (1973)  (Ref.  51  developed  a  new  method  of 
calculating  both  the  complex  propagation  constant  t  =  jk  and  characteristic 
impedance  of  a  material  using  transmission  line  theory  applied  to  an  impedance 
tube  device.  His  method  consists  of  measuring  standing  wave  properties  in  the 
proximity  of  an  absorbing  material  alternately  backed  with  a  high  impedance 
and  a  low  impedance  termination. 

Finally,  in  1980,  Chunq  and  Blaser  [Ref.  61  published  their  article 
discussing  a  Dual-Microphone  Transfer  Function  technique  for  measuring  in-duct 
acoustic  properties,  namely  the  acoustic  impendance  and  reflection 
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coefficient  from  which  the  complex  wave  number  can  be  determined.  This 
technique  was  compared  to  the  conventional  Standinq  Wave  Tube  technique  and 
found  to  qive  similar  results  [Ref.  7}. 

C.  PURPOSE  OF  STUDY 

The  objective  of  this  research  is  to  investigate  the  feasibility  of 
obtaininq  the  characteristic  impedance  and  the  complex  wave  number  of  a  fluid- 
filled  porous  solid  by  measuring  the  resonance  frequencies  and  half-power 
widths  of  a  fluid-filled  rectangular  cavity  which  is  partially  occupied  by  a 
porous  solid  of  varying  thickness.  In  this  investigation,  the  fluid  is  air, 
the  porous  solid  is  fiberglass  and  the  cavity  is  made  of  plexiglass. 

This  investigation  differs  from  previous  investigations  in  that  1)  Biot 
Theory  [Ref.  8,  9,  10,  and  11]  is  applied  to  solve  for  the  allowed  modes  of 
the  partially-filled  cavity,  and  2)  for  the  primary  mode  of  interest  the 
particle  velocity  is  directed  along  the  length  of  the  fluid-porous  solid 
boundary  as  opposed  to  normal  to  this  boundary. 

Table  1  surma rizes  the  advantages  in  the  course  of  this  research.  The 
advantage  of  this  technique  over  the  currently  existing  methods  is  that  it 
does  not  require  probing  the  sound  field  either  within  or  in  front  of  the 
porous  solid  as  do  methods  1,2,  3,  and  5,  nor  does  it  require  a  pair  of 
phase-matched  microphones  as  does  method  6.  Additionally,  characteristic 
properties  of  an  absorbing  naterial  need  not  be  known  as  in  method  4.  Also, 
the  acoustic  properties  of  the  bulk  fiberglass  can  be  determined  using  only 
thin  samples. 

The  theory  developed  for  the  new  method  is  presented  in  Chapter  II  along 
with  a  general  overview  of  the  equations  used  in  the  classical  methods.  The 
experimental  results  for  each  of  the  methods  used  in  comparison  are  presented 
in  Chapter  III  and  comparisons  are  made  and  discussed  in  Chapter  IV. 
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TABLE  1 


ADVANTAGES  AND  DISADVANTAGES  OF  METHODS  USED  10  DETERMINE  THE  CHARACTERISTIC 

IMPEDANCE  AND  COMPLEX  WAVE  NUMBER 


Method 

1.  Probe  Tube 
(Scott) 


2.  Impedance  Tube 
( Beranek ) 


Advantages 
a)  Simple  theory. 


b)  Simple  procedure, 
a)  Simple  theory. 


b)  Simple  procedure. 

c)  Probing  of  material 
not  required. 

d)  Normal  absorption  easy 
to  measure. 


Disadvantages 

a)  Thick  sample  needed 
to  get  character¬ 
istic  impedance 
without  a  wave 
model  for  the 
porous  medium. 


b)  Must  probe 
f ield . 


sound 


a)  Thick  sample  needed 
to  get  character¬ 
istic  impedance 
without  a  wave 
model  for  porous 
medium. 


3.  Resonance 

a) 

Probing  not  required. 

a) 

Theory  is 

somewhat 

(Beranek) 

difficult. 

b) 

Measurement  of  re¬ 

b) 

Discrete 

frequen- 

sonance  frequency  and 
half  power  widths  is 

cies. 

precise. 

c) 

Measures 

impedance. 

surface 

4.  Material  Properties 
( Beranek ) 


d)  Requires  low  losses 
for  well  defined 
resonances . 

a)  Sound  wave  not  needed.  a)  Complex  theory. 

b)  Measurement  of 
microscopic  proper¬ 
ties  is  complex. 
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TABLE  1  (continued) 


Method 


Advantages 


Disadvantages 


c)  Requires  classifi¬ 
cation  of  material 
as  rigid  or  soft. 


5.  Yaniv 


a)  Simple  theory. 


a)  Tine  consuming. 


b)  Simple  procedure. 


b)  Quarter  wavelenqth 
termination  needs 
continuous  adjust¬ 
ment  with  frequen¬ 
cy. 


6.  IVo  Point  Phase  Dif- 
erence  (Chung  and 
B laser) 


c)  Probing  of  material 
not  necessary. 

a)  Simple  theory. 


b)  Simple  procedure. 


a)  Thick 
needed. 


sample 


b)  Microphone  posi¬ 
tions  only  optimum 
for  select  freguen- 


c)  Probinq  of  material 
not  necessary. 


c)  Requires  accurate 
microphone  nhase 
calibrations. 


7,  Mew  Method 
( Baker) 


a)  Probinq  not  required. 


b)  Simple  procedure. 


a)  Complex  theory  and 
numerical  methods 
required. 

b)  Discrete  frequen¬ 
cies. 


c)  Ability  to  determine 
DC  flow  resistance 
non-invasively. 


c)  Requires  low  losses 
for  well  defined 
resonances . 


d)  Bulk  properties 

obtained  from  thin 
samples. 


II .  THEORY 


The  theory  behind  the  new  method  is  presented  in  detail  in  this  section 
with  much  of  the  fundamental  concents  and  understanding  havinq  come  from  the 
investiqation  of  sound  nropaqation  in  a  porous  solid  saturated  with 
super-fluid  Helium  [Ref  11].  Additionally,  the  basic  theory  and  equations  used 
are  presented  for  each  of  the  classical  methods  used  in  comparison  with  the 
new  techniques. 


A.  THEORY  OF  A  FLUID-SATURATED  RIGID  POROUS  SOLID 

The  linearized  equation  of  motion  tor  the  microscopic  fluid  velocity 
vnicro  the  Eulerian  description  is  [Ref.  12,  p.  49], 


Pf 


3v, 


micro 


3t 


+  ^Vmicro*  ^vmicro 


=  -5?  +  n72vnicro  +  U  +  ^  n)V(V-v 


micro 


(2.1) 

). 


Using 


+■  +  -*• 

y2V  .  =  V(  v.V  .  )  _ 

v  micro  vv  micro' 


*  *  * 

7  y  V  v  V  . 

v  *  v  *  micro 


(2.2) 


we  can  write  Eq.  (2.1)  as 


Pf 


i^icro  +  (4ioco.;,;ntcrct 


*  -7n  +  iW  x  7  x  vmicro  +  ( C  +  3  n)7(  ?*vnicro^ ' 


Where  pf  is  the  fluid  mass  density,  p  is  the  fluid  pressure,  n  is  the  fluid 
shear  viscosity  and  ?  is  the  fluid  bulk  viscocity.  The  last  term  in  this 
equation  can  be  ingnored  here  because  its  affect  is  small  compared  to  the 
other  terms  on  the  scale  of  a  pore.  Ignoring  the  last  term,  Eq.  (2.3)  becomes 


Pf 


Dv, 


micro 


Dt 


=  -Vp  -  n^kVxv, 


micro 


(2.4) 


Dv^icro 

where  — -  is  the  total  time  deriviative  of  the  microscopic  fluid 

velocity.  Fran  Eq.  (2.4),  a  phenomenological,  linearized  equation  for  the 
average  (over  the  microscopic  structure)  fluid  velocity,  denoted  <vmicro>'  03,1 
be  written: 

peff  -  ^<P>  -  Kflow<Vcro>  •  (2.5) 

Here  is  the  phenomenological  effective  fluid  mass  density  and  Rjclow  is 
the  phenomenological  flow  resistance  (Ref.  13,  pp.  252-255).  Rfiow  is  defined 
such  that  in  steady  state 

-7<p>  *  Rdc  <*micro>  .  (2.6) 

In  general  (including  oscillatory  flow)  pgff  and  Rfj^  can  be  operationally 
defined  by 


1  <V 

2  Peff  v  micro' 

*  hK?51Sof3h!gy  per 

(2.7) 

♦  2 

Rflow  <Vmicro> 

*  Rate  of  dissipation  of 
kinetic  energy  density 

(2.8) 

For  ej“*-  time  dependance,  Eq.  (2.5)  can  be  written,  upon  combining  both 
velocity  terms  on  the  left  hand  side, 

%ff(“)  +  Rf low( <Vmicro>  *  '  *  (2*9) 
where  ogff  and  Rfiow  are  now  frequency  dependent  phenomenological  properties 
which  describe  the  effect  on  the  fluid  due  to  the  microscopic  geometry  (Ref. 
II,  pp.  334-335].  It  should  be  noted  that  in  Eq.  (2.9)  and  in  the  remainder  of 
this  thesis  complex  quantities  are  denoted  by  a  tilde. 
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Defininq  complex  effective  density,  peff(SJ),  as 

'■*  **  R f-  ]  ow (  w) 

Peff(u)  =  Peff{u)  + - -  )  (2.10) 

Eq.  (2.9)  may  be  written  in  a  more  compact  form  as 

jwPeff^)  <Vnicro>  =  ~7<p>  .  (2.11) 

Plane  proqressive  wave  solutions  to  Eq.  (2.11)  travellinq  in  the  positive 
x-direction  nay  be  souqht  by  lettinq  the  acoustic  pressure  be  qiven  by 
5p  “  kx^  and  lettinq  < ''micro'**  ^avq  wt  ”  kx^  where  k  =  k  -  ja  with 

k,  a  >  0  such  that  the  enerqy  of  the  wave  exponentially  decays  as  one 
moves  deeper  into  the  porous  solid  (note  in  order  to  avoid  confusion  later, 
the  letter  a  will  be  used  to  denote  the  attenuation  constant,  reserving 
the  letter  a  for  the  tortuosity  in  the  Biot  theory).  Eq.  (2.11)  then 
becomes 


mm 


where  cf  is  the  speed  of  sound  in  the  bulk  fluid,  ^U)  is  the  complex  ratio 
of  specific  heats  (Cj/cv)  as  a  function  of  frequency,  and  y  is  the  adiabatic 
ratio  of  specific  heats. 

Eq  s.  (2.12),  (2.13)  and  (2.14)  lead  to  a  wave  equation  with  the  phase 
velocity  determined  as  follows:  writing  Eq.  (2.12)  (the  equation  of  notion) 


Jf  Peff<w)  vavq  =  6p 


and  substituting  Eq.  (2.14)  into  it  results  in 


(2.15) 


f  Peff<^  vavg  =  <=f2  «p  . 


(2.16) 


Recasting  Eq.  (2.13)  (continuity)  as  pg  v 
by  it  results  in 


Pf  *avg  “  £  5p  and  dividing  Eq.  (2.16) 


y 

y 

i* 

v 

f  I 

» 


I 

C' 


(2.17) 


Defining  the  complex  tortuosity  in  the  fluid  saturated  porous  solid  by  (Ref. 
8,  p.  185  and  Ref.  11,  p.  37] 


A»ff(B> 


the  phase  velocity  can  be  written  as 


(2.18) 


£2  *  (Z)2  s  f£ _  i(il 

The  ocmplex  tortuosity  can  be  written  as 


a(  w)  *  a  + 


RqqF( w) 


(2.19) 


(2.20) 
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The  quantity  r/ct  determines  the  relative  importance  of  flow  resistance  versus 
fluid  inertia  in  detemininq  the  microscopic  flow  field. 

The  structural  S  depends  upon  the  averaqe  cross-sectional  shape  of  a 
pore  [Ref.  8,  p.  185]  .  Experimentally  it  appears  that  6=4  for  many  qranular 
porous  solids  [Ref.  16] . 

Substitutinq  Eq.  (2.22),  (2.23)  and  (2.24)  into  Eq.  (2.21)  results  in 

(t)2  (i§r)]1/2l  •  U.26) 

For  all  calculations  in  this  research  it  will  be  assumed  that  a=l  and  5=4  for 
fiberqlass  [Ref.  11,  pp.  49-59]. 

The  complex  ratio  of  specific  heats,  y(  u>)  ,  is  defined  as 

y(uj)  =  Cp/cv  (2.27) 

where  cp  and  cv  are,  respectively,  the  complex,  frequency  dependent  specific 
heats  per  qram  at  constant  pressure  and  at  constant  volume.  For  an  ideal  qas 

^  ^  R 

Cp  -  Cy  =  —  ,  where  M  is  the  molar  mass,  and  so 


*")  =1  +  M5v 


(2.28) 


The  propaqation  of  sound  in  the  porous  medium  will  vary  from  isothermal 
at  low  frequencies  when  the  thermal  penetration  depth  is  larqe  compared  to 
a  pore  size,  hence  cv  +  *•  ,  to  adiabatic  at  hiqh  frequencies  when  the  thermal 
penetration  depth  is  small  compared  to  a  pore  size,  and  cv  =  bulk  fluid 
value.  Since  the  transport  of  momentum  and  enerqy  in  a  qas  both  arise  from 
molecular  collisions,  it  will  be  assumed  that  thefunctional  form  (y  -  1)  has 

<v  <v 

the  same  form  as  the  complex  tortuosity  a(w),  only  with  the  thermal 
penetration  depth  5^h  in  place  of  the  viscous  penetration  depth  5visc. 
Hence  we  take  [S.  Baker,  private  cormunication] 


yU)  =  1  +  1-j  *1(3) 


(2.29) 


Sx 
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where  r  and  F(  w)  are  given  by  Eq's  (2.22)  and  (2.23),  respectively.  The 
viscous  penetration  depth,  <$visc  '  ^  the  thermal  penetration  depth  <5^  are 
defined  by 


=  /HT  5th  =  (2,30) 

where  n  is  the  fluid  shear  viscosity,  pf  is  the  bulk  fluid  density,  <  is 
the  thermal  conductivity  and  Cp  is  the  specific  heat  per  gram  at  constant 
pressure. 

6th 

The  factor  — — -  can  be  recast  as 
'Vise 


where 


(2.31) 


and 


v  *  Vpf  =  kinematic  viscosity 


(2.32) 


x  3  </p  Cp  =  thermal  diffusivity  (2.33) 
and  Pr  =  Prandtl  number  which  has  a  value  of  0.733  for  air  at  20°C.  Hence  the 
complex  ratio  of  specific  heats  can  be  rewritten  as 

■ 1  +{^rfF(a  •  <2-341 

TT 

In  the  high  frequency  limit,  the  square  of  the  phase  velocity  in  the 
porous  material  (Eq.  2.19)  qoes  to  the  fluid  inertia  dominated  limit  in  which 
sound  propagates  at  the  bulk  fluid  speed  reduced  by  the  square  root  of  the 
tortuosity.  Conversely,  in  the  low  frequency  limit,  flow  resistance 
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dominates,  and  pressure  oscillations  diffuse  throuqh  the  porous  material.  The 
crossover  from  inertia  dominance  to  flow  resistance  dominance  occurs  when  the 
viscous  penetration  depth  is  approximately  equal  to  a  characteristic  pore 
size.  Figure  2.1  sumnarizes  the  limiting  conditions. 

1 .  Development  of  Dispersion  Relation  in  a  Fluid-Filled  Cavity  Partially 


Filled  With  a  Porous  Solid 


Consider  a  rectangular  cavity  which  is  partially  filled  with  a  fluid- 
saturated  porous  solid.  The  region  occupied  by  bulk  fluid  is  denoted  as 
region  1;  the  region  occupied  by  the  fluid-saturated  porous  solid  is 
identified  as  region  2.  The  porous  solid  is  assumed  homogeneous  and 
isotropic.  Figure  2.2  depicts  the  geometry  used  to  develop  the  dispersion 
relation  for  standing  waves  in  the  cavity.  It  will  be  assumed  that  all  cavity 
cavity  walls  are  perfectly  rigid  so  that  n*u  *  0.  Then  n»Vp  =  0  at  these 
boundaries  and  the  velocity  potential  in  regions  1  and  2  can  be  written  as 
follows; 

In  region  1  (bulk  fluid). 


►  1  =  <(>io  e',U)t  cos(kxx]  cos[klz(z-d1)] 


(2.35) 


where 


nir 


k  =  —  ;  n  =  0,  1,  2,  .  .  .  , 

X  L 


w2  *  (kx2  +  klz2)  Cf' 


(2.36) 

(2.37) 


.2  = 


Bf/pf 


(2.38) 


and 


Bf  =  bulk  modulus  of  fluid 

$10  =  complex  amplitude  of  the  velocity  potential  in  region  1. 


23 


VJ 

vj 


8 

$ 


VJ 

Vj 

5 


! 


£ 


•TJ 


3 


.V 


21 

r.] 

t 

A’ 


a. 


’  i 

V 

yi 

1 


ww^varawaaiMBffltii^^ 


High  and  Low  Frequency  Limits  on  Phase  Velocity 
in  a  Porous  Material 


High  Frequency  Limit  Low  Frequency  Limit 


Plug  Flow  Laninar  Flow 


o(  2)  +  au» 

y(u»)  ♦  y 


al  ai)  *  - 

Y(  u)  ♦  1 


Jpf3 


Cf' 


.  Cf2upf 


Fluid  Inertia  Dominates 
Sound  Propagates 


Flow  Resistance  Dominates 

Sound  does  not  Propagate 
(Pure  Diffusion) 


Cross-over  when: 

5yisc  -  characteristic  pore  size 


Figure  2.1  High  and  Low  Frequency  Limits  on  Phase  Velocity  in  a  Porous 
Material 


REQXON  2.  C^OROUS  MATtRlA^ 


Figure  2.2  Rectangular  Cavity  Geometry  for  New  Method 


K> 


The  x-conponent  of  the  wave  nunber  is  a  pure  real  nunber  since  the  boundary 
conditions  of  riqid  walls  normal  to  the  x-direction  at  x  =0  and  x  =  L  force 
real  and  discrete  values  for  kx  . 

In  reqion  1  the  particle  velocity  can  be  written  as 


vi  -  7*i 


(2.39) 


and  the  pressure  as 

Spl  -  ~ju  Pf  ?i 

In  reqion  2  (porous  material), 

A-  =  A _ d  1  AAC  f  U 


where 


k  = 


nn 


x  L 


and 


(2.40) 


cos  [  k  22  (  z+d  2 )  ] 

(2.41) 

,2 ,  .  .  .  , 

(2.42) 

.  2  ~  2 

■  2Z  , 

(2.43) 

y(  w) 

w  » 

(2.44) 

<{>20  =  complex  amplitude  of  the  velocity  potential  in  reqion  2. 
The  particle  velocity  in  reqion  2  is 


V  2  ~  v  4>  2 


(2.45) 


and  the  pressure  is,  from  Eq.  (2.11), 


With  the  velocity  potentials  defined  in  both  regions  1  and  2  all  that 
is  left  is  to  natch  conditions  at  the  boundary  between  the  bulk  fluid  and 
porous  material  (i.e.,  at  z  *  0)  to  produce  the  desired  dispersion  relation. 
The  two  boundary  conditions  that  must  be  satisfied  at  z  =  0  are: 

1.  continuity  of  pressure  across  the  boundary: 

4g(z=o)  =  a( Z)  $2( z=o)  (2.47) 

and 

2.  continuity  of  the  normal  component  of  volume  velocity: 


34>1 

Tt~ 

2*0 


P 


3  $2 

TT" 

z=o 


(2.48) 


where  P  is  the  porosity  of  the  material.  Apolyinq  boundary  condition  1 
results  in 


$10  costk^dj  *  a(u>)  $20  cos [k 2zd2] 


(2.49) 


Applying  boundary  condition  2  results  in 

$io  *iz  sin  [k izdt}  *  -$2o  P  k^  sin(k2zd2]  .  (2.50) 

Dividing  Eg.  (2.50)  by  Eg.  (2.49)  results  in  the  dispersion  relation 


where 


klz  tan  [k12d!)  +  k2z  tan  [k2zd2]  =  0  ,  (2.51) 


•iz 


U» 


=£2  - 


1/2 


[2.52) 
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V-  U  tA  '< 


<***  a;*'*** 


2 


1/2 


(2.53) 


and 


^  “<“>  W  kx 

f 


n  =  0,  1,  2,  .  .  . 


(2.54) 


In  Eq.  (2.53)  the  forms  of  a(w)  and  y(w)  have  already  been  defined  in  Eq. 
(2.26)  and  Eq.  (2.34),  respectively.  Eq.  (2.51)  is  a  transcendental  equation 
in  a)  whose  roots  determine  the  natural  frequencies  and  half  power  widths  or, 
equivalently,  the  quality  factors,  of  the  allowed  "standing  wave"  modes  of 
this  system.  Fran  these  roots  can  be  determined  the  complex  wave  numbers 

i'v 

and  k2z  and  the  resulting  pressure  distribution. 

Two  approaches  to  the  use  of  the  dispersion  relation  become  evident. 
The  first  is  the  "forward"  problem,  where  the  roots  of  Eq.  (2.51)  for  u>  may 
be  regarded  as  functions  of  Rqc  anc*  a  ,  and  measured  resonance  frequencies 
and  0's  may  be  compared  with  predictions;  the  second  in  the  "inverse" 
problem,  where  Rqc  and  <*  are  regarded  as  functions  of  the  roots  of  Eq. 
(2.51)  for  w  and  are  found  by  variation  such  that  the  best  fit  is  obtained 
between  the  measured  and  calculated  spectra.  The  aim  of  this  work  is  to 
verify  the  technique  to  solve  the  inverse  problem  by  comparing  the  value  of 
RoC  extracted  from  resonance  data  with  that  actually  measured.  Appendix 
A  discusses  the  numerical  technique  to  determine  the  roots  of  the  dispersion 
relation  and  the  algorithm  used  to  determine  the  best  fit  value  of  Rpc* 

Once  the  "best  fit"  value  of  Rq£  has  been  determined,  then,  since 
the  roots  of  the  dispersion  relation  in  terms  of  u>  are  known  for  this  value 
of  Roc  ,  the  complex  wave  number  is  found  from 
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(2.56) 


- - 


with  c  found  frcm  Ecj's.  (2.19),  (2.26)  and  (2.34).  The  characteristic 


impedance  can  be  calculated  from 


^char  =  feff ^  c 


(2.57] 


where  /eff(to)  is  given  by  Eq.  (2.10).  Noramlizing  the  characteristic 
impedance  to  that  of  the  bulk  fluid  (^qc)  results  in 


Zchar  (normalized)  =  Zchar//»0c 


(2.58) 


The  complex  pressure  field  for  both  regions  may  be  calculated  once 


the  roots  of  Eq.  (2.51)  are  found.  The  pressure  distribution  in  region  1  can 


be  written  in  terms  of  the  velocity  potential  in  region  1  as 


£p,  =  -P/’f  ?i  ' 


(2.59) 


where  <f>  ^  is  given  by  Eq.  (2.35).  The  pressure  distribution  in  region  2  can 


be  written  in  terms  of  the  velocity  potential  in  region  2  as 


S p,  ■  “j55  *(S)/f  , 


(2.60) 


where  ^2  *-s  given  by  Eq.  (2.41). 

The  velocity  potentials  and  02  are  related  by  the  boundary 


condition  that  pressure  be  continuous  across  the  interface  between  the  bulk 


fluid  and  the  porous  material,  Eq.  (2.47),  which  requires  that 


Ao  cos[k,  zd,  1  =  «(u>)  Via  cos[k2zdz]  . 


(2.61) 


For  simplicity,  choose  ^20  to  depend  on  jt>10  >  resulting  in 


2  -  cos  k|Zd| 

20  “  *(«5)  cos  1<2Z  dz  * 


(2.62) 


The  pressure  distribution  in  each  region  may  be  written  as  follows: 


.  ** 
-Dw 


ft  cos(k  x]  cosfk,  z(z-d, )] 


I  K 


(2.63) 


.V  .VA  J>  A  A  A.VA.V.V^ 


and  for  region  2, 


LCOS  K  .  Q  _ 

-ju  Pf  cos[kxx]  —  s'^2zd^  cos[k2z(z+d2)]  ?10  eJ<iJ  (2.64) 

*>* 

With  Eq.  (2.63)  and  Eq.  (2.64),  and  the  values  of  kj2  and  k2Z  determined 
from  the  roots  of  the  dispersion  relation  of  Eq.  (2.51),  the  pressure 
distribution  throughout  the  cavity  may  be  calculated.  See  Appendix  B  for  a 
more  detailed  discussion  of  the  pressure  calculation. 

B.  CURRENT  METHODS  USING  TRANSMISSION  LINE  THEORY 

Three  of  the  five  classical  methods  for  determining  the  characteristic 
impedance  and  complex  wave  number  of  a  porous  acoustic  material  are  based  on 
transmission  line  theory.  These  methods  involve  measurement  of  various 
properties  of  the  standing  wave  and  data  are  taken  at  various  points  somewhat 
removed  from  the  sample  surface  (Ref.  4,  pp.  317-3211.  The  three  methods 
considered  ares 

—  Impedance  Tube  Method 
—  Resonance  Analysis  Method 
—  Yaniv's  Method 

Only  Yaniv's  method  allows  direct  determination  of  the  material's 
characteristic  impedance  and  complex  wave  number.  The  impedance  tube  method 
allows  direct  measurement  of  the  resistive  and  reactive  components  of  the 
surface  impedance.  The  oomplex  wave  number  can  only  be  calculated  if 
reflections  from  the  far  side  of  the  porous  material  can  be  neglected  and  one 
has  some  knowledge  of  the  complex  effective  mass  density  of  the  fluid  in  the 
pores.  As  developed  by  Beranek,  the  resonance  anlaysis  method  measures  the 
same  quantities  as  the  impedance  tube  method,  with  the  same  limitations.  It  is 


discussed  primarily  for  its  similarity  to  the  new  technique  in  that  the 
parameters  of  the  resonance  curve  are  measured. 

1 .  Impedance  Tube  Method 

In  this  method  the  non-dissipative  acoustic  transmission  line  is 
considered  a  smooth  rigid-walled  tube  whose  lateral  dimensions  are  small 
compared  to  a  wavelength.  The  losses  at  the  wall  and  in  the  bulk  volume  are 
considered  snail  compared  to  the  losses  associated  with  the  absorbing  material 
under  test.  Assuming  that  the  source  of  sound  is  at  one  end  of  the  tube  (x=0) 
and  the  sample  under  test  is  located  at  the  opposite  end  (x=L),  the  surface 
impedance  can  be  determined.  Figure  2.3  depicts  the  geometry  of  the  impedance 
tube. 

The  impedance  tube  that  has  a  constant  cross  section  S  and  the 
source  vibrates  harmonically  at  a  frequency  sufficiently  law  that  only  plane 
waves  can  propagate.  A  standing  wave  pattern  is  set  up  by  the  right  and  left 
travelling  waves  and  can  be  expressed  mathematically  as 

p  ,  A  ejIut  +  tBe’1'*  -  ka-x)l  (2.65) 

where  p  is  the  acoustic  pressure,  A  and  B  are  determined  from  the  boundary 
conditions  at  x=0  and  x=L. 

It  has  been  shown  [Ref.  17,  pp.  200-206]  that  the  surface  impedance 
of  the  termination  may  be  written  in  the  form, 
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(2.66) 


where  B/A  is  the  ratio  of  the  relected  wave  amplitude  to  the  incident  wave 
amplitude  and  9  is  the  phase  difference  between  the  relected  and  incident 
wave.  Thus  given  B/A  and  9  ,  the  surface  impedance  of  the  termination  can  be 
calculated. 
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Solving  for  the  amplitude  |p|  of  the  standing  wave  it  can  be  seen  that 
the  amplitude  at  pressure  antinode  is  A+B,  and  at  a  pressude  node  is  A-B.  The 
ratio  of  the  pressure  amplitude  at  an  antinode  to  that  at  a  node  is  the 


standing  wave  ratio,  or  SWR: 


SVJR  = 


A+B  ^  Panti 

Pnode  ’ 


which  can  be  rearranged  to  give  the  ratio  B/A: 

B  SWR-1 
A  *  SWR+1  * 


(2.67) 


(2.68) 


Thus  measurement  of  the  SWR  by  probing  the  sound  field  in  the  impedance  tube 
with  a  snail  microphone  yields  the  value  of  B/A  (Ref.  17,  p.  206] . 
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Figure  2.3  Geometry  for  the  Impedance  Tube 


The  phase  angle  9  can  be  evaluated  from  determining  the  distance  to 


the  first  node  from  the  end  at  x*L.  These  nodes  are  located  at 
k(L-xn)  -  9/2  ■  (n-1/2) w  ;  n  ■  1,  2,  .  .  .  , 

so  that  the  phase  angle  for  the  first  node  is 

0i  *  2k(L-x ,)  -  it  . 


(2.69) 


(2.70) 


With  the  SWR  and  8^  known,  the  resistive  and  reactive  portions  of  the 
impedance  can  be  determined  graphically  using  a  "Smith"  chart.  To  determine 
the  characteristic  impedance  of  the  porous  material,  the  sample  thickness  is 
increased  such  that  the  reflected  wave  from  the  rigid  termination  is 
insignificant  compared  to  the  incident  wave.  From  the  definition  of  the 
material's  characteristic  impedance,  the  complex  wave  number  can  be  found 
assuming  that  we  know  how  the  fluid  reacts  in  the  pores  of  the  material.  The 
characteristic  impedance  may  be  written  as 

2 char  s  peff  c  (2.71) 

where  pgff  is  the  complex  effective  fluid  density  in  the  pores  and  c  is  the 

phase  velocity  of  the  acoustic  disturbance  in  the  pores.  The  phase  speed 

<s, 

c  in  a  fluid,  can  expressed  as 

c2  *  Bf/Pef f  (2.72) 

where  ^  is  the  complex  bulk  modulus  of  the  fluid  in  the  pores.  Therefore 
Eq.  (2.71)  becomes,  upon  substituting  Eq.  (2.72)  for  Peff' 

>v  -v 

Zchar  *  Bf/c  (2.73) 

and  since  c  =  the  equation  that  describes  the  complex  wave  number  It  can 
be  written  as 

k  *  Zchar  •  (2.74) 

If  is  is  assumed  that  the  material  pore  size  is  small  compared  to  the  viscous 
and  thermal  penetration  depths  (in  air  these  two  parameters  are  approximately 
equal)  in  the  frequency  range  of  interest,  we  can  write 

Bf  -  &T  (2.75) 

where  Bp  is  the  isothermal  bulk  modulus  of  the  fluid  in  the  pores. 
Additionally  for  an  ideal  gas,  the  isothermal  bulk  modulus  can  be  written  as 
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&T  =  l/ioj.  (2.76) 
where  is  the  isothermal  compressibility  of  the  gas,  and  is  equal  to  1/P, 
where  P  is  the  pressure  in  the  pores  [Ref.  18,  p.  47).  Hence  Eq.  (2.74)  can 
be  written,  assuming  isothermal  conditions  exist,  in  the  form, 


k  —  ujiop  Zchar  (2.77) 
where  k  *  k  -  ja,  for  a  wave  travelling  in  the  positive  x-direction  in  the 
porous  solid.  Eq.  (2.77)  is  presented  to  show  the  relationship  between  the 
complex  wave  number  and  the  characteristic  impedance  assuming  isothermal 
conditions  in  the  pore  volume.  This  equation  was  used  only  to  check  Beranek's 


Material  Properties  Method  in  the  isothermal  limit.  Only  the  values  for 
’Zcftar  extracted  by  the  impedance  tube  method  are  reported  here. 

2.  Resonance  Analysis  Method 

The  resonance  analysis  method  is  used  to  determine  the  surface 
impedance  of  a  porous  acoustic  material.  This  method  requires  low  loss  within 
the  material  to  be  successful.  A  resonance  curve  may  be  traced  out  by  varying 
either  the  tube  length  or  the  frequency  in  this  method.  For  the  present 
research,  the  impedance  tube  discussed  previously  was  used,  with  frequency 
swept  to  trace  out  a  resonance . 

Beranek  presents  the  equation  needed  to  determine  the  impedance  of  a 
sample  under  test  and  is  rewritten  here  only  for  discussion  [Ref.  4,  pp.  329- 


334] . 


The  normalized  surface  impedance  of  the  sample  is  written, 


lsu£f  3  ¥/*  -i  kj/i^  ,  coth  -  jo^t/c] 

PqC 


(2.78) 


where  ¥  =  |ZSurfl/poc  '  $  being  the  phase  angle  of  |Zsurf|  and  Zsurf 
is  the  impedance  at  the  surface  of  the  sample  which  is  to  be  determined.  At 
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the  resonance  peak,  the  peak  frequency  ojp  is  approximately  equal  to  the 
normal  angular  frequency  of  longitudinal  vibration,  ^  ,  or 


“p  5  4n  • 


12.79) 


At  the  half  power  frequencies  (-3  dB  down  from  the  peak)  the  half-power  half 
width  (HPHW)  can  be  written  as 


HPHW  =  —= - -  =  — — -  =  kl  , 

z 


(2.80) 


where  kj  =  +  k^,  and  kna  is  the  spatial  damping  constant  without  the 
porous  material  present  (i.e.,  it  takes  into  account  wall  and  driver  losses). 
Therefore  1^,  the  spatial  damping  factor  due  to  the  sample  under  test  can  be 


written  as 


k  =  k,  _  k 


(2.81) 


Assuming  the  ratio  kn/%  «  1,  the  surface  impedance  of  the  sampole  can  be 
determined  once  u^,  u^,  %,  and  c  are  measured. 

3.  Yaniv’s  Method 

Equations  have  been  developed  in  transmission  line  theory  for 
determination  of  the  input  impedance  of  an  electrical  transmission  line  in 
terms  of  the  line  characteristics:  propagation  constant,  characteristic 
impedance  and  length  of  line  [Ref.  51 .  When  the  line  is  terminated  in  an 
open  or  short  circuit,  the  input  impedance  of  the  line  may  be  expressed  in 
terms  of  the  line  parameters  with  the  following  equations: 


2 os  =  zchar  ^^h  tyLl 


(2.82) 


zss  =  zchar  f  Y^l 


(2.83) 


where  2^  is  the  input  impedance  to  a  line  terminated  in  an  open  circuit 
(Figure  2.4),  Zss  is  the  input  impedance  to  a  line  terminated  in  a  short 


35 


mmrnaaaamaatmaanamtBrvi  mjh  *r*  rar 


wngganuncn  »n  mr\  «--i  «n  wi  Kn 


circuit  (Figure  2.5),  ZChar  *s  the  characteristic  impedance  of  the  line, 
Y  =  jk  is  the  complex  propagation  constant  of  the  line  and  L  is  the  length 
of  the  line. 

From  Eq.  (2.82)  and  Eq.  (2.83)  we  can  solve  for  the  material's 
characteristic  impedance: 

Zchar  -  Uss  W  1/2  '  (2.34) 

and  for  the  complex  propagation  constant 

tanh(YL]  *  tZgg/Z^og]  1/2  .  (2.85) 

The  high- impedance  termination  for  an  acoustical  line  may  be  realized  by 
termination  with  a  rigid  wall  and  a  low-impedance  termination  for  an 
acoustical  line  may  be  realized  by  terminating  the  line  with  a  column  of  fluid 
equal  to  a  quarter  wavelength,  itself  terminated  in  a  rigid  wall  [Ref.  51 . 


Figure  2.4  Geometry  for  Yaniv's  Open  Circuit  Determination 


Figure  2.5  Geometry  for  Yaniv’s  Short  Circuit  Determination 
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By  measuring  the  standinq  wave  inside  the  tube  alternatively 
terminated  with  an  open  circuit  and  then  a  short  circuit,  both  zos  and  zsg 
may  be  determined.  With  these  two  values,  Eq.  (2.84)  and  Eq.  (2.85)  may  be 
solved  for  the  desired  complex  propagation  constant  and  characteristic 
impedance. 

C.  PROBE  TUBE  METHOD 

This  method  was  first  presented  by  Scott  (1946)  [Ref.  2]  to  allow  the 
determination  of  the  complec  wave  number  from  the  characteristics  of  an 
acoustic  disturbance  propagating  through  a  fluid-saturated  porous  solid  by 
probing  the  sound  field  existing  in  the  material  and  extracting  the  wave 
number  and  absorption  coefficient. 

To  determine  the  complex  wave  number  in  the  porous  medium,  which  describes 
the  spatial  variation  of  the  acoustic  disturbance  for  a  uniform,  progressive 
plane  wave,  Scott  writes  the  velocity  potential  $  of  the  average  motion  in  the 
pores  as, 

♦  »  e^6*  (2-86) 
where  k  =  k  -  ja  is  the  complex  wave  number  and  x  represents  distance  in  the 
direction  of  propagation  of  the  wave.  Eq.  (2.86)  may  be  rearranged  such 
that, 

~  =  toe"35'  "  kx)  (2.87) 
where  a  is  the  absorpotion  coefficient  and  k  *  2tt/X  is  the  wave  number. 
Experimentally,  one  measures  the  acoustic  pressure  anplitude  as  a  function  of 
distance  in  the  material,  along  with  the  wave  length. 

To  determine  the  characteristic  impedance  of  the  porous  material,  Scott 


writes, 


where  Peff^)  is  the  complex  frequency  dependant  effective  density  of  the 
fluid  in  the  pores.  Determination  of  Zchar  assumes  of  course,  that  the 
effective  fluid  density  is  known,  for  example,  from  Biot  theory.  Eq.  (2.88) 
was  used  only  to  check  that  the  measured  values  of  k  were  consistent  with  the 

'■v 

values  of  ZChar  from  the  other  methods.  Only  the  values  of  k  measured  with 
Scott's  method  are  reported  here. 

D.  MATERIAL  PROPERTIES  METHOD 

This  method  requires  determination  of  five  properties  of  the  material  to 
calculate  the  complex  propagation  constant  and  characteristic  impedance  of  a 
material.  Beranek  [Ref.  4,  Ch.  19]  discusses  this  method  in  detail  and  only 
that  portion  of  the  theory  applicable  to  the  comparison  with  the  new  technique 
will  be  presented. 

Reranek  has  shown  that  both  the  propagation  constant  and  characteriustic 
impedance  are  derivable  with  the  aid  of  charts  from  the  basic  physical 
properties  of  the  material  [Ref.  4,  p.  839]. 

These  properties  are: 

(a)  the  specific  flow  resistance 

(b)  the  porosity 

(c)  the  structure  factor 

(d)  the  volume  coefficient  of  elasticity  of  the  air  in  the  interstices 

(e)  the  volime  coefficient  of  elasticity  of  the  skeleton  of  the  material. 
First,  one  must  determine  whether  the  material  under  investigation  can  be 
considered  a  soft  blanket  or  a  rigid/dense  blanket.  This  determination  is 
made  by  considering  the  ratio  of  the  volume  coefficient  of  elasticity  of  the 


air  in  the  interstices  to  that  of  the  skeleton.  If  this  ratio  is  less  than 
Yq  the  material  is  considered  to  be  a  riqid/dense  blanket  [Ref.  4,  p.  840] . 
The  fiberglass  material  chosen  for  the  rigid/dense  blanket  criterion  for 

<v 

Beranek  provides  an  equation  for  the  propagation  constant  b  =  jk,  of  such  a 
material,  namely 


b  =  jk  =  jw 


'PQkst^ 


P0kstw 


where 

b  =  complex  propagation  constant  (same  as  Yaniv's  y) , 
k  =  complex  wave  number, 

kst  =  structure  factor  (same  as  Biot's  tortuosity,  a), 

K  =  volume  coefficient  of  elasticity  of  air  in  the  interstices  (i.e., 
the  bulk  modulus,  and  can  be  isothermal,  adiabatic  or  complex), 

Rj  =  dynamic  specific  resistance  per  unit  thickness  (approximately  the  DC 
flow  resistance), 

Pq  =  denisity  of  air,  and 

Y  =  porosity  (same  as  P  in  Biot's  theory). 

By  the  use  of  Figure  2.6  and  Figure  2.7  the  magnitude  and  phase  of  b  can  be 

determined.  The  characteristic  impedance  may  then  be  determined  from 


2  .  =  ~j  ^77  b 

char  J  wY 


[2.90) 


This  method  suffers  from  the  fact  that  not  all  of  the  material  properties  can 
be  experimentally  determined  and  thus  reasonable  assumptions  for  their  values 


must  be  made  [Ref.  4,  p.  840] . 
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III.  EXPERIMENT  AND  RESULTS 


The  objective  of  this  study  was  to  develop  a  new  technique  to  treasure  both 
the  complex  wave  number  and  characteristic  impedance  of  a  fluid-saturated 
porous  solid.  As  with  any  new  technique,  comparison  with  classical  methods  is 
a  qood  basis  on  which  to  judge  one's  theory.  This  section  will  first  present 
the  experimental  results  of  the  five  classical  techniques  discussed  previously, 
and  then  the  new  technique.  A  oonparison  of  measured  data  for  all  methods  will 
be  presented  in  Chapter  IV. 

A.  PROPERTIES  OF  THE  POROUS  90LID 

The  porous  solid  used  in  all  experiments  was  a  semi-rigid  type  of 
fiberglass,  much  like  Owens-Coming  Tupe  705  fiberglass,  used  in  the  insulation 
of  refrigeration  type  equipment.  This  material  has  sufficient  rigidity  to 
allow  it  to  be  cut  from  sheets  into  a  variety  of  shapes  and  to  allow  it  to  be 
inserted  into  tubes  without  excessive  compression.  This  property  was  important 
in  the  flow  resistance  and  impedance  tube  measurements  where  a  snug  fit  was 
required. 

Discs  with  nominal  dimensions  3-29/32  inches  in  diameter  and  1-29/32 
inches  thick  were  cut  from  three  sheets  of  this  material,  weighed,  and  measured 
to  determine  the  material's  volume.  From  two  sets  of  measurements,  the  average 
material  density  was  found  to  be  43.0  t  5.0  kg/m3.  For  comparison, 
Owens-Coming  Type  705  fiberglass  is  reported  to  have  a  denisty  of  96  kg/m3 
(Ref.  5),  while  rock-wool,  also  known  commercially  as  Stillite,  is  reported  to 
have  a  density  of  80  kg/m3  [Ref.  2]. 

The  porosity  of  the  fiberglass  material  was  measured  by  iflinersinq  a  1-inch 
cube  in  a  graduated  beaker  of  water  and  measureing  the  volume  of  displaced 
water.  The  fiberglass  material  itself  was  assumed  to  be  non-absorbinq.  The 

resulting  porosity  was  found  to  lie  between  0.99  and  1.0. 
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1.  Flow  Resistance  Measurements 


Next,  measurement  was  made  of  the  specific  DC  flow  resistance  per  unit 
lenqth  of  sample.  For  most  homogeneous,  isotropic  materials  this  flow 
resistance  may  be  measured  with  sufficient  accuracy  by  applyinq  a  known  steady 
air  flow  across  a  sample  and  measuring  the  resulting  pressure  differential. 
The  actual  method  used  is  similiar  to  that  presented  by  Beranek  [Ref.  4,  p. 
850] . 

The  quantities  measured  were:  [Ref.  4,  p.  847] 
ip  =  pressure  differential  between  the  two  faces  of  the  sample 
v  =  velocity  of  the  linear  flow  of  air  through  the  sample  due  to  &p 
The  pressure  differential  was  measured  with  a  standard  water  manometer  (±0.1 
cm  accuracy)  and  the  velocity  of  linear  flow  was  determined  buy  measuring  the 
time  required  to  collect  a  sufficient  volume  of  water  (which  drew  air  through 
the  sample  by  a  siphon)  in  a  graduated  cylinder.  The  experimental  apparatus 
used  is  shown  in  Figure  3.1.  All  tube  connections  were  sealed  with  RTV  Silicon 
Rubber  to  prevent  leaks. 

The  velocity  of  linear  flow  produced  by  ip  is, 


where 


(3.1) 


V  =  volume  of  water  collected  in  time  t,  and 
A  =  cross-sectional  area  of  sample. 

Then  the  DC  flow  resistance  Rq^  is  given  by 

Rrv,  =  M  =  A£At  =  P^AhAt 
«DC  vT  VT  VT 


(3.2) 


where  p  is  the  density  of  water,  g  is  the  acceleration  of  gravity,  Ah  is 

the  difference  in  water  height  in  the  two  arms  of  the  manometer,  and  T  is  the 

thickness  of  the  sample.  The  units  of  Rqc  ace  in  Ray 1/unit  thickness,  which 

can  by  written  in  GGS  or  MKS  units  gnV’cm2,sec/an  or  kg/m2,sec/m,  respectively. 
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In  the  first  flow  resistance  experiment  seven  discs  were  placed  into 
the  larqe  sample  holder,  each  disc  havinq  a  diameter  of  9.85  an  and  a  thickness 
of  4.76  an.  The  samples  were  made  as  close-fittinq  as  possible  but, 
unfortunately,  the  die  used  to  cut  these  discs  had  a  diamter  approximately  0.1 
cm  greater  than  that  of  the  sample  holder.  While  providing  a  tight  fit  along 
the  tune  wall,  this  led  to  a  compression  in  the  material's  thickness  by 
approximately  10%  as  they  were  pushed  into  the  holder.  At  that  time  it  was 
thought  that  the  effect  of  this  compression  would  not  be  significant  enough  to 
require  re-machining  the  cuting  die.  Later  measurements  of  flow  resistance 
using  a  different  holder,  however,  gave  values  which  were  a  factor  of  2  snaller 
than  that  for  flow  normal  to  the  fiberglass  surface. 

Measurements  were  oonducted  over  a  limited  range  of  air  flow  velocity 
by  adjusting  the  3/4  inch  Hoke  valve  on  the  water  side  of  the  apparatus  to  give 
the  desired  pressure  differential.  Mien  steady  state  conditions  were  reached, 
the  volume  flow  was  then  plotted  as  a  function  of  volume  flow  rate.  This  is 
shown  in  Figure  3.2.  Within  experimental  error,  the  variation  in  Ap  with  V  was 
found  to  be  linear  up  to  the  highest  flow  rate  measured,  which  was  many  orders 
of  magnitude  greater  than  expected  for  a  sound  wave.  Since  the  sound  pressures 
appearing  at  the  surface  of  a  sample  are  seldom  greater  than  1000  dynes/an2  in 
practice  (100  Pa  or  0.5  inches  water)  it  can  be  seen  that  even  at  higher 
pressure  differenitals  the  linear  relationship  betwen  Ap  and  V  is  maintained. 
The  DC  flow  resistance  had  a  value  of  59,100  ±  4,200  MCS  Ray 1/m  using  the 
uncompressed  thickness.  Table  2  lists  the  experimental  data  and 
Table  3  lists  the  DC  flew  resistance  results  using  the  compressed  and 
uncompressed  thickness. 

It  was  discovered  later  (while  attempting  to  reconcile  these 
measurements  with  the  DC  flow  resistance  required  to  fit  the  resonance  data) 
that  the  DC  flow  resistance  of  our  fiberglass  samples  is  quite  anisotropic  (the 


TABLE  2 


LARGE  SAMPLE  HOLDER  DC  FLOW  RESISTANCE  P1EASUREMENT  DATA 


(Flow  Normal  to  Surface) 


Uncompressed  thickness  =  33.34  ±0.1  cm 
Compressed  thickness  =  29.80  ±0.1  cm 


Sample  holder  diameter  *  9.84  ±0.1  cm 
Rocm  tanperature  =  23.3  deg.  Centriqrade 


Ah  Vol. 

(cc)  (cc) 


6.4  ±  0.1 

2038 

6.4 

2102 

6.3 

2074 

5.6 

2016 

5.5 

2056 

5.3 

2008 

4.5 

2046 

4.5 

2122 

4.3 

2108 

3.5 

2128 

3.5 

2063 

3.5 

2074 

2.5 

2077 

2.5 

2062 

2.5 

2103 

1.7 

2098 

1.7 

2081 

1.7 

2072 

1.1 

2086 

1.2 

2120 

1.2 

2100 

0.6 

2110 

0.6 

2116 

0.6 

2161 

t 

(sec) 

7.3  ±  0.1 

7.6 

7.4 

8.4 

8.7 

8.9 

10.3 
10.6 

11.4 

14.1 

14.3 

13.5 

18.1 
18.8 
19.0 
27.0 
26.8 

26.7 

46.9 
46.0 

44.3 

98.5 

94.4 
93.3 


V 

(cc/sec) 

278.0  ±  0.5 

277.7 

278.8 

240.6 

236.6 

225.1 

198.6 

201.1 

184.9 

151.5 

144.7 

153.6 

114.6 

109.5 

110.9 

77.9 

77.6 

77.5 

44.4 

46.1 

47.4 

21.4 

22.4 

23.2 


Roc* 

(MKS  Ray 1/m) 

57.5x10 3  t  1 

57.6 

56.5 

58.2 

58.1 

58.8 

56.6 

55.9 
60.8 

57.7 

60.4 

56.9 

54.5 

57.1 

56.3 

54.5 

54.8 

54.8 

61.8 

65.1 

63.3 
70.0 
66.8 

64.7 


*Ccmpressed  Thickness  Data 


TABLE  3 


LARGE  SAMPLE  HOLDER  DC  FLOW  RESISTANCE  MEASUREflENT  RESULTS 

(flow  normal  to  surface) 

Compressed  Thickness; 

Mean  =  59.1  x  10 3  MKS  Ray  1/m 
Std.  Deviation  =  4.2  x  10 3  MKS  Rayl/m 

Uncompressed  Thickness; 

Mean  Rqq  =  53.2  x  10 3  MKS  Rayl/m 
Std.  Deviation  =  4.2  x  10 3  MKS  Rayl/m 
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^  smaller  sanr>le  holder  was  subsequently  constructed  and  samples  were  cut  which 
allowed  the  measurement  of  flow  resistance  with  particle  velocity  nomal  to  the 
fiberqlass  surface  and  with  particle  velocuty  parallel  to  the  fiberqlass 
surface.  Figure  3.3  shows  the  results. 

Mot  only  is  the  DC  flow  resistance  anisotropic  by  a  factor  of  two  to 
three  but  the  effect  of  the  10%  compression  of  the  material  in  the  larqe  sample 
holder  is  evident.  The  resultinq  flow  resistance  parallel  to  the  surface  of 
the  material  was  found  to  be  15,200  t  500  MKS  Ray 1/m  while  that  normal  to  the 
surface  was  34,000  t  900  MKS  Rayl/m.  Table  4  and  Table  5  qive  the  measured 
data. 


2.  Viscous-Inertial  Transition  Frequencv 


The  ratio  of  the  viscous  penetration  depth  to  the  characteristic  pore 


size  is  one  of  the  most  fundamental  dimensionless  parameters  that  determines 
the  nature  of  elastic  wave  propoqation  in  a  fluid-saturated  porous  solid.  If 
Vise  <<:  P°re  size»  then  the  notion  of  the  fluid  in  the  pores  is  essentially 
the  sane  as  if  the  fluid  were  inviscid.  However,  if  5vjsc  >>  pore  size  then  the 
fluid  in  the  pores  is  essentially  frozen  to  the  solid,  and  sound  cannot 
propaqate  within  the  fluid.  Microscopic  photoqraphs  of  the  fiberqlass  material 
provided  an  estimate  of  the  averaqe  distance  between  fibers.  This  mean 
distance  was  determined  to  be  220  ±  80 in.  The  crossover  from  where  fluid 
inertia  dominates  to  where  flow  resistance  dominates  occurs  when  6visc  «*  1/2 
(mean  pore  diameter).  This  value  was  estimated  to  be  110  t40  gm.  Since,  the 
thermal  penetration  depth  in  air  is  approximately  the  same  as  the  viscous 
penetration  denth,  the  crossover  from  isothermal  to  adiabatic  compression  in 
the  pores  was  expected  to  be  nearly  coincident  with  the  viscous- inertial 
transition.  Figure  3.4  shows  the  viscous  penetration  depth  over  the 
frequencies  of  interest.  The  transition  frequency  for  this  case  was 
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approximate!- '  100  :'z.  for  a  collection  of  straight  oioes  of  radius  a,  tue  x; 
flow  resistance  is  given  ov, 

Rnc  =  (3.3' 

Considerin'!  the  T  flow  resitance  measured  wits  flow  parallel  to  tne 
air- fiberglass  surface,  the  eguivalent  radius  would  oe  9-!  ^m.  This  fiiuo- 

compares  very  well  with  the  nicroscooic  observations.  \  radius  of  )H  _m  r°sjlts 
in  a  transition  freouency  of  500  Hz. 


h.  CLASSIC  'KTHODS  OF  TRANSMISSION  LINE  IMEORy 
1 .  Impedance  Tube  Experiment  and  Results 

The  impedance  tube  technigue  is  the  most  widely  used  to  determine  a 
material's  characteristic  impedance,  being  relatively  simple  in  noth  theory 
and  measurement.  This  experiment  was  the  first  conducted  of  the  five  classical 
methods. 

A  Rruel  and  Kjaer  Type  4002  Standing  Wave  apparatus  was  used  for  these 
measurements.  The  tube  was  constructed  of  clear  Plexiglass  which  allowed 
visual  inspection  of  the  sample.  The  tube  is  1  meter  in  length  and  9.0  cm  in 


diameter.  This 

sets  a 

lower  limit  on 

the 

usable  freouency  range  to 

be 

approximately  35 

Hz  for 

the  length  of 

the 

tube  to  be  greater  than 

1/4 

wavelength,  and  a  higher  limit  of  approximately  2006  Hz  for  operation  below  the 

lowest  crossmode,  thus  ensuring  only  plane  waves  are  propagated.  the  HP4192A 

Low  Freouency  Impedance  Analyzer  was  used  for  signal  generation  and  detection. 

It  has  a  built  in  freouency  synthesizer  with  a  resolution  of  1  mHz  and  a 

phase-sensitive  detector  with  a  0,1  Hz  bandwidth  to  preclude  noise  and 

harmonics  interfering  with  measurements.  Figure  3.5  shows  the  experimental 

setup.  Since  we  are  primarily  interested  in  the  lowest  longitudnal  mode  in  the 

rectangular  cavity  (particle  velocity  parallel  to  the  surface  of  the 

fiberglass),  which  has  a  freouency  in  the  range  of  600  to  700  Hz,  the  impedance 
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Press.  Difference  vs  Volume  Flow 

(FLOW  NORHRL/'PRRflLLEL  TO  SURFRCE) 


VOLUME  FLOW  RATE 
(ee/i*e ) 


Figure  3.3  Ap  vs  V  in  the  Small  Sample  Holder 


Rvg .  Pore  Size  &.  Viscous  Depth 

(TRANSITION  REGION  LIMITS) 


FREQUENCY 

(Hi) 


Figure  3.4  Characteristic  Length  Scale  of  a  Pore  Compared  to  the  Viscous 
Penetration  Depth 
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TABLE  4 


SMALL  SAMPLE  HOLDER  DC  FLOW  RESISTANCE  MEASUREMENT  DATA  AND  RESULTS 

(Flow  Normal  to  Surface) 

Thickness  =  15.60  ±0.01  cm 
Cross-Sectional  Area  =  20.5  ±0.1  cm2 


Ah 

Vol. 

t 

V 

Rdc 

(cc) 

(cc) 

(sec) 

( cc/sec ) 

(MRS  Rayl/m) 

7.6  ±  0.1 

2165  ±  4 

7.8  ±  0.1 

278.6  ±  0.5 

35.1x10 3  ±  900 

7.5 

2086 

7.6 

275.2 

35.1 

6.0 

2052 

9.0 

229.3 

33.7 

5.9 

2096 

9.4 

222.3 

34.2 

4.9 

2082 

11.2 

185.6 

34.0 

5.0 

2108 

10.9 

192.9 

33.4 

4.0 

2084 

13.5 

154.3 

33.4 

4.0 

2116 

13.4 

157.4 

32.7 

3.0 

2150 

18.3 

117.6 

32.9 

3.0 

2098 

18.0 

116.6 

33.2 

2.0 

2118 

28.4 

74.6 

34.6 

2.0 

2151 

29.1 

74.0 

34.8 

1.0 

2116 

57.8 

36.6 

35.2 

1.0 

2121 

55.4 

38.3 

33.7 

Mean  Rqq  =  34.0  x  10  3  MRS  Ray  1/m 
Std  Deviation  =  0.9  x  10 3  MRS  Rayl/m 


TABLE  5 


SMALL  SAMPLE  HOLDER  DC  FLOW  RESISTANCE  ‘1EASUREMENT  CATA  AND  RESULTS 

(Flow  Parallel  to  Surface) 

Thickness  =  17. BO  ±  0.01  on 
Cross-Sectional  Area  =  20.5  t  0.1  cn2 


Ah 

Vol. 

t 

V 

Rdc 

(cc) 

(cc) 

( sec) 

( cc/sec ) 

(MKS  Rayl/n) 

3.6  ±  0.1 

2068  t  4 

7.5  ±  0.1 

276.8  ±  0.5 

14.7x10 3  t  500 

3.6 

2100 

7.5 

281.1 

14.4 

3.1 

2066 

9.0 

229.8 

15.2 

3.1 

2120 

9.1 

232.2 

15.1 

2.5 

2109 

11.9 

177.7 

15.9 

2.5 

2124 

11.1 

191.5 

14.7 

1.9 

2126 

14.9 

142.8 

15.0 

1.9 

2115 

15.2 

139.0 

15.4 

1.3 

2075 

22.1 

94.0 

15.6 

1.3 

2108 

22.9 

92.0 

15.9 

0.7 

2150 

42.2 

51.0 

15.5 

0.7 

2104 

39.6 

53.3 

14.9 

Mean  =  15.2  x  10 3  MKS  Rayl/m 

Std.  Deviation  =  0.5  x  10  3  flKS  Rayl/n 
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tube  measurements  were  carried  out  at  frequencies  between  100-1000  Hz  with 
sample  thicknesses  of  4,  8,  and  12  inches.  The  different  sample  thicknesses 
were  used  in  order  to  determine  how  thick  the  material  had  to  be  for 
reflections  from  the  rigid  termination  to  be  negligible,  thus  allowing  the 
material's  characteristic  impedance  to  be  determined. 

Measurements  of  the  reflection  coefficient  from  100  to  1000  Hz  for  all 
three  samples  indicated  that  the  material  could  be  considered  infinite  for 
thicknesses  greater  than  8  to  12  inches.  The  reflection  coefficient,  R,  is 
given  by 

R  =  i  B/A.  |  ,  (3.4) 
where  the  ratio  8/A  is  given  by  Eq.  (2.68).  The  reflection  coefficient  chanqes 
noticeably  in  the  frequency  region  betwen  400  and  500  Hz.  This  is  thouqht  to 
be  due  to  the  transition  from  isothermal  to  adiabatic  compression  of  the  gas  in 
the  pores. 

The  standing  wave  ratio  and  the  distance  to  the  first  pressure  node 
were  measured  for  the  12"  sample  at  frequencies  between  100  and  1000  Hz.  The 
normalized  characteristic  impedance  !zChar!/poc  was  calculated  usinq  the 
"Smith"  chart  nomoqram  with  the  SiJR  and  phase  at  the  first  node  as  input 
parameters.  Table  6  gives  the  results. 

It  was  observed  that  at  frequencies  below  200  Hz  fluctuations  occur  in 
both  magnitude  and  phase  of  the  impedance.  The  position  and  the  amplitude  of 
the  pressure  antinode  and  node  were  very  difficult  to  measure  accurately  below 
200  Hz,  and  so  these  data  are  not  as  reliable  as  the  data  taken  above  200  Hz. 
From  200  to  400  Hz,  the  impedance  magnitude  decreases  at  a  larqer  rate  than  at 
frequencies  above  500  Hz,  apparent  evidence  of  the  isothermal  to  adiabatic 
transition  occurring  in  the  pores.  This  transition  is  seen  in  both  magnitude 
and  phase  as  well  as  the  reflection  coefficient.  The  transition  coincides 
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Figure  3.5  Experimental  Setup  for  the  Impedance  Tube  Method 
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TABLE  6 


IMPEDANCE  TUBE  MEASUREMENT  RESULTS  ON  12”  OF  FIBERGLASS 


Freq. 

(Hz) 

SWR 

Distance 
to  first 
node (cm) 

A# 

1 ^char 1/ Roc 

$7 

£char 
( Deqrees ) 

100 

2.79 

59.70 

1.60 

-44.0 

125 

5.83 

59.70 

3.74 

-48.0 

200 

4.79 

36.80 

2.64 

-52.7 

250 

4.40 

28.75 

2.87 

-45.6 

300 

4.00 

24.10 

2.83 

-42.1 

400 

3.46 

18.00 

2.60 

-38.0 

500 

2.23 

14.45 

1.95 

-22.9 

600 

2.47 

11.40 

1.94 

-31.7 

700 

2.26 

9.65 

1.80 

-29.3 

800 

2.30 

8.25 

1.75 

-31.7 

900 

2.15 

7.75 

1.81 

-24.4 

1000 

2.11 

6.90 

1.77 

-25.1 
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with  the  velocity  profile  chanqinq  from  laminar  flow  in  the  low  frequency  limit 
to  pluq  flow  in  the  hiqh  frequency  limit  (see  Fiqure  2.1). 

2.  Resonance  Analysis  Experiment  and  Results 

Of  all  the  classic  mehtods,  this  technique  in  most  similar  to  the  new 
technique  in  that  parameters  are  extracted  from  the  resonance  curve  to 
determine  a  sample’s  characteristic  impedance.  As  previously  cemented,  this 
method  is  limited  to  discrete  frequencies  and  thin  samples  since  the  resonances 
of  the  tube  become  weaker  as  thicker  samples  are  introduced.  The  upper  limit 
thickness  for  our  fiberqlass  was  found  to  be  about  2  inches.  A  1  inch  thick 
sample  was  chosen  since  Yaniv's  method  was  conducted  with  the  same  sample 
thickness  and  since  three  of  the  five  empty  tube  resonances  were  still  distinct 
enouqh  to  make  measurements. 

The  experimental  setup  is  very  similar  to  the  impedance  tube 
experiment.  A  Bruel  and  Kjaer  Type  4002  Standinq  Wave  apparatus  was  used  alonq 
with  a  HP3580A  Spectrum  Analyzer.  The  spectrum  analyzer  was  used  to  determine 
the  resonance  frequencies  and  half-power  widths  of  the  lowest  three  peaks  in 
both  the  empty  tube  and  in  the  tube  filled  with  the  1  inch  thick  fiberqlass 
sanple.  A  HP3314A  Function  Generator  was  then  set  to  the  frequency  of  each  of 
these  peak  and  the  distance  between  the  fiberqlass  surface  and  the  pressure 
maximum  closest  to  the  speaker  was  measured.  This  is  needed  to  calculate  the 
effective  tube  length,  as  an  infinite  impedance  loudspeaker  is  assumed  in  the 
theory.  The  surface  impedance  of  the  sample  was  then  determined  from  these 
measurements  usinq  Eq.  (2.78).  Fiqure  3.6  shows  the  experimental  setup  and 
Table  7  gives  the  resultinq  values  of  normalized  impedance  values. 

The  lowest  two  resonances  (center  frequencies  of  253  Hz  and  400  Hz, 
respectively)  qave  plausible  values  of  surface  impedance,  however  the  highest 
mode  with  the  absorbing  material  present  (center  frequency  of  538  -  544  Hz)  had 
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TABLE  7 


RESONANCE  ANALYSIS  SURFACE  IMPEDANCE  RESULTS 
(Impedance  Tube  and  Rectanqular  Cavity) 

(1  Inch  Thick  Fiberglass  Sample) 


Surface  Impedance  Results  in  Impedance  Tube 


Resonance 

Peak 

Number 

Center 

Frequency 

(damped) 

(Hz) 

HPFW 

(Hz) 

1  ^ surf  1 / Poc 

£ 

zsurf 
( degrees ) 

1 

253 

23 

2.9 

-85.9 

0.0078 

2 

400 

18 

2.1 

-88.9 

0.0013 

3* 

543 

41 

7.7 

+66.4 

-0.0055 

Surface  Impedance  Results  in  Rectangular  Cavity 


Resonance 

Peak 

Number 

Center 

Frequency 

(damped) 

(Hz) 

HPFW 

(Hz) 

*0 

1 ^surf 1 / Poc 

^surf 
( deqrees ) 

iSn^Sa 

1 

841 

126 

1.8 

-65.4 

0.071 

*  This  point  was  discarded,  as  discussed  in  text. 
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a  higher  0  than  in  the  empty  tube  and  was  therefore  discarded.  The  erroneous 
data  for  the  third  resonance  peak  is  attributed  to  tne  significant  broadening 
of  its  resonance  curve  and  the  resultant  overlapping  with  adjacent  nodes. 

This  resonance  analysis  method  was  applied  to  the  rectangular  cavity 
with  and  without  a  1  inch  thick  strip  of  fiberqlass  while  exciting  the  [010] 
mode  to  verify  the  results  obtained  in  the  tube.  The  [010]  mode  was  chosen 
because  for  this  mode  the  particle  velocity  is  normal  to  the  fiberglass 
surface,  as  in  the  impedance  tube.  For  these  measurements  in  the  rectangular 
cavity  the  source  was  placed  at  grid  point  (0,0)  (Figure  3.11)  and  the  receiver 
at  grid  point  (20,0).  The  resonance  curve  was  first  measured  in  the  empty 
cavity;  it  had  a  center  frequency  of  873  Hz.  The  normalized  magnitude  of  the 
surface  impedance  obtained  agreed  closely  with  the  value  obtained  by  the 
impedance  tube  for  frequencies  between  800  and  900  Hz.  However  a  significant 
phase  discrepency  was  evident.  The  reason  for  the  discrepency  is  not  known. 
Perhaps  it  is  a  result  of  using  a  point  source  located  in  a  corner  .  Table  7 
lists  the  results. 

As  with  other  impedance  tube  measurements,  this  method  relies  upon  a 
plane  propagating  wave,  where  incidence  on  the  porous  material  is  assumed 
normal  to  the  surface  of  the  sample  under  test.  With  a  point  source  located  in 
one  corner  of  the  cavity  used,  departure  frcrn  plane  wave  propagation  can  be 
expected  with  increasing  sample  thickness  as  the  quality  factor  of  the 
resonance  degrades.  The  extracted  impedances  are  surface  values.  In  both 
geometries  in  which  this  method  was  used  the  magnitude  of  the  surface  impedance 
decreased  with  increasing  frequency  (which  was  expected).  The  phase  in  all 
measurements  indicates  that  the  impedance  remains  largely  reactive,  and  that 
the  qas  in  the  pores  acts  like  a  stiffness.  Use  of  samples  of  qreater 
thickness  would  only  suppress  the  resonance  curves  more. 
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.  Yaniv's  ‘tethod,  Experiment  and  Results 

The  third  and  last  technique  which  uses  transmission  line  theory  is  the 
method  proposed  and  experimentally  verified  by  Yaniv  [Ref.  5] .  The  experiment 
was  conducted  usinq  the  impedance  tube  of  the  previous  two  experiments,  with 
several  nodif ications.  Yaniv's  method  requires  neasurinq  the  surface  impedance 
with  the  sample  backed  by  a  riqid  termination  and  by  an  acoustic  short-circuit, 
obtained  by  movinq  the  riqid  termination  such  that  a  1/4  X  air  qap  exists 
between  the  rear  surface  of  the  material  and  the  riqid  termination. 

The  riqid  end  cap  used  in  the  two  previous  experiments  was  replaced 
with  a  2  inch  thick  disc  of  plexiqlass  with  a  3/8  inch  diameter  hollow  aluminum 
tube  threaded  into  the  side  away  from  the  driver.  A.  small  amount  of  hiqh 
vacuum  qrease  was  used  to  seal  the  edges  of  the  disc  at  the  tube  walls  while 
still  allowinq  for  positioning.  The  front  face  of  the  disc  had  a  recess  cut 
out  of  it  into  which  a  FET  microphone  mounted  at  the  end  of  a  3/16  inch  hollow 
aluminum  rod  could  fit  snuqly  and  flush  so  as  to  provide  a  gnooth  surface  when 
the  disc  was  placed  aqainst  the  back  end  of  the  sample.  The  3/16  inch  diameter 
rod  fit  inside  the  3/3  inch  diameter  rod,  providinq  support  and  allowinq  the 
microphone  to  be  positioned  anywhere  between  the  backside  of  the  porous  solid 
and  the  riqid  termination.  Electrical  connection  to  a  power  supply  switchbox 
was  also  made  by  wires  which  passed  throuqh  the  small  hollow  rod.  This 
conf iquration  allowed  the  disc  and  microphone  to  serve  as  the  riqid  termination 
for  the  open  circuit  impedance  measurement  and  also  allowed  accurate 
positioning  of  the  termination  to  provide  the  required  one-quarter  wavelength 
air  gap  between  the  porous  solid  and  the  riqid  termination  at  each  frequency. 
Fiqure  3.7  shows  the  movable  termination  and  Figure  3.8  shows  the  equipment  set 
setup. 

A  1  inch  thick  disc  of  fiberglass  was  inserted  into  the  impedance  tube 
such  that  the  distance  betwen  the  driver  and  the  sample  was  at  least  a  quarter 


Figure  3.7  The  Rigid  Sliding  Termination  (Yaniv'S  Method) 

wavelength.  The  lowest  usable  frequency  was  therefore  approximately  250  Hz, 
the  upper  limit  remaining  approximately  2000  Hz  to  preclude  the  lowest 
crossmode.  Measurements  were  made  of  the  first  minimum  and  maximum  pressure 
amplitudes  in  front  of  the  sample  and  of  the  distance  between  the  sample 
surface  and  the  first  pressure  minimum  in  from  of  the  sample.  For  each 
frequency,  these  measurements  were  performed  under  conditions  of  open  and  short 
circuit  terminations.  A  short  circuit  termination  at  the  back  side  of  the 
fiberglass  disc  was  obtained  by  adjusting  the  position  of  the  plexiglass  plug 
until  a  null  output  was  observed  from  the  microphone,  which  was  held  flush  to 
the  fiberglass  disc.  It  was  found  that,  within  experimental  error,  the  result¬ 
ing  distance  between  the  porous  disc  and  the  plexiglass  plug  was  within  1% 
of  the  1/4  wavelength  determination  made  in  the  empty  tube  at  the  frequencies 
used  in  the  experiment.  The  standing  wave  ratio  and  the  reflection  coefficient 
phase  angle  were  calculated  using  Eq.  (2.67)  and  Eg.  (2.70),  respectively.  The 
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open  and  short  circuit  impedances  were  found  using  the  "Smith"  chart.  The 
characteristic  impedance  was  then  found  by  use  of  Eq.  (2.84). 

Determination  of  the  complex  propagation  constant,  y  =  jk,  required 
finding  the  roots  of  Eq.  (2.85)  such  that, 

tanhfyL]  -  [Z^g/Z^I  1  /2  =  0  .  (3.5) 
The  roots  were  found  using  a  numerical  technique  known  as  Ward's  method.  The 
details  are  given  in  Appendix  A  [Ref.  19] .  The  results  of  this  experiment  are 
listed  in  Table  8. 

As  for  the  previous  methods,  a  decrease  in  impedance  magnitude  with 
increasing  frequency  was  observed,  with  the  transition  from  isothermal  to 
adiabatic  compression  between  400  to  700  Hz  in  this  1  inch  thick  sample.  The 
complex  wave  number  also  shows  this  transition.  The  data  and  results  for  400 
Hz  appear  to  anomalous.  The  reason  is  not  known. 

C.  PROBE  TUBE  METHOD,  EXPERIMENT  AND  RESULTS 

The  probe  tube  measuremnets  were  made  by  setting  up  a  plane  wave  in  the 
Bruel  and  Kjaer  Type  4002  Standing  Wave  apparatus  filled  partially  with 
fiberglass  and  exploring  the  sound  dield  by  means  of  a  probe  tube  extension  to 
the  impedance  tube  microphone.  The  slope  of  the  decay  of  the  pressure 
amplitude  with  distance  provided  a  measurement  of  the  attenuation  constant  a. 
The  change  in  phase  of  the  oscillatory  wave  with  distance  provided  the  wave 
number  k.  The  oomplex  wave  number  is  then  found  from  k  =  k  -  ja. 

The  equipment  setup  is  shown  in  Figure  3.9.  Three  specimens,  of 
thickness  4,  8,  and  12  inches,  were  chosen  for  measurement  and  comparison  with 
the  impedance  tube  results.  Each  sample  thickness  was  made  up  of  the  required 
number  of  fiberglass  discs,  9.84  cm  diameter  and  4.76  cm  thick,  inserted  into 
the  tube  opposite  the  driver,  taking  care  to  avoid  compression  of  each  disc 
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TABLE  8 


YANIV'S  METHOD  EXPERIMENTAL  RESULTS 


2os 


Jss 


? 


Freq. 

(Hz) 

SVJR 

Distance  to 
first  node 
(cm) 

SWR 

Distance  to 
first  node 
(cm) 

250 

50.40 

30.90 

1.36 

30.50 

300 

42.40 

25.80 

1.33 

25.70 

400 

24.03 

18.60 

1.89 

33.20 

500 

24.71 

14.10 

1.68 

21.70 

600 

12.29 

11.50 

1.65 

17.20 

700 

9.57 

9.70 

1.56 

14.70 

800 

8.71 

8.20 

1.70 

12.60 

900 

7.16 

6.80 

1.70 

11.15 

1000 

6.42 

5.90 

1.67 

10.45 

Freq. 

(Hz) 

lzcharl/P0c 

^zchar 
( deqrees ) 

k 

(Rad/m) 

a 

(Np/m) 

cph 

(m/s) 

250 

2.79 

-44.8 

12.82 

13.96 

122.5 

300 

2.78 

-43.6 

12.56 

13.86 

150.2 

400 

2.04 

-22.3 

15.17 

8.80 

165.7 

500 

2.17 

-30.1 

20.82 

12.75 

150.9 

600 

2.10 

-28.9 

21.70 

15.77 

173.7 

700 

1.98 

-28.4 

21.73 

16.81 

202.2 

800 

1.94 

-27.9 

24.88 

17.50 

202.1 

900 

1.75 

-27.7 

27.75 

18.40 

203.7 

1000 

1.62 

-25.6 

28.41 

17.92 

221.2 

&&<»">  y?  Xh/i toy?*'  jCt  /.v.ylv  / 
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while  trying  to  minimize  any  air  gaps  between  disc  surfaces.  Each  disc  had  a 
0.71  cm  diameter  hole  cut  in  the  center  for  a  close  fit  around  the  probe  tube 
without  any  deformation  of  the  fiberqlass. 

The  probe  tube  was  manufactured  out  of  steel  tubing  with  the  same  outside 
and  inside  diameter  of  the  existing  cpen  ended  probe  tube  on  the  standing  wave 
apparatus.  A  copper  cap  was  welded  onto  the  end  of  the  extension  next  to  a 
3/16  inch  diameter  hole  to  allow  the  measurement  of  the  sound  field.  The 
extension  attached  to  the  existing  probe  tube  with  a  plastic  sleeve.  This  is 
shown  in  Figure  3.10. 

The  experiment  was  conducted  at  frequencies  of  125,  250,  500,  and  1000  Hz 
for  each  of  the  three  sample  thicknesses,  but  only  the  12  inch  thick  sample 
was  long  enough  so  that  the  wave  reflected  from  the  rigid  end  of  the  tube  was 
sufficiently  attenuated  to  be  negligible  compared  to  the  incident  wave.  The 
pressure  amplitude  at  the  rigid  termination  was  measured  to  be  20  dB  less  than 
that  at  the  sample  boundary.  Determination  of  the  attenuation  constant  was 
made  by  plotting  pressure  versus  distance  on  semi-log  paper  and  extracting  the 
slope.  The  phase  of  the  pressure  was  measured  with  a  SR-350  Dual  Channel 
Lock-In  Amplifier.  The  wave  number  was  determined  from  the  change  of  phase 
with  distance  into  the  fiberglass.  Table  9  displays  the  experimental  results 
in  terms  of  k  and  a. 

The  impedance  magnitude  decreases  with  frequency  while  the  phase  indicates 
that  the  material  becomes  more  resistive  with  increasing  frequency.  Only  the 
absorption  coefficient  shows  any  sign  of  the  transition  between  isothermal  and 
adiabatic  compressions  in  the  pore  volume. 

D.  MATERIAL  PROPERTY  METHOD,  EXPERIMENT  AND  RESULTS 

All  the  methods  discussed  so  far  involve  measuring  certain  characteristics 
of  a  progressive  plane  wave  in  the  porous  material  or  its  proximity.  Beranek 
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Figure  3.9  Equipment  Setup  for  Probe  Tube  Experiment 


has  shown  that  the  complex  wave  number  and  characteristic  impedance  can  also  be 
derived  from  the  basic  physical  properties  of  the  material  with  the  aid  of 
charts.  These  properties  are: 

a)  specific  flow  resistance  per  unit  length,  Rj(a  Rqq) 

b)  porosity,  Y 

c)  structure  factor,  kst  (tortuosity) 

d)  volume  coefficient  of  elasticity  of  air  in  the  interstices  K  (Bulk 
modulus) 

e)  volume  coefficient  of  elasticity  of  the  skeleton  of  the  material. 
Measurements  of  the  specific  flow  resistance  and  porosity  have  already  been 
discussed.  The  remaining  three  material  properties  are  difficult  if  not 
impossible  to  measure,  and  will,  based  on  a  reasonable  knowledge  of  the 
material,  be  approximated. 

The  structure  factor  kst,  known  as  the  tortuosity  in  Biot  theory,  is 
essentially  a  path  length  correction  factor.  Experimentally,  it  is  found  to 
depend  most  strongly  on  the  porosity  of  the  porous  solid  [Ref.  11] .  In  the 
limit  porosity  *  1,  the  tortuosity  also  equals  1  for  typical  porous  maerials, 
and  certainly  for  fiberglass.  As  the  porosity  of  the  samples  used  in  this 
research  was  found  to  be  very  close  to  1,  a  value  of  1  for  the  tortuosity  was 
assumed  in  all  calculations. 

The  volume  coefficient  of  elasticity  of  air  in  the  interstices,  K,  is  an 
important  term;  if  isothermal  conditions  can  be  assumed  then  K  2  1  x  10 6  dynes/ 
cm2;  if  adiabatic  conditions  exist  then  K  2  1.4  x  10 6  dynes/cm2.  For  most 
fibrous  acoustic  materials  however,  K  will,  in  same  frequency  reqion,  be 
neither  isothermal  nor  adiabatic  but  will  be  seme thing  in  between  and  will  be 
complex  valued  and  frequency  dependent,  as  is  our  case. 

Finally,  the  volume  coefficient  of  elasticity  of  the  skeleton,  Kskel' 
although  not  directly  involved  in  the  calculation  of  k  and  Zchar,  must  be 
known  or  approximated  since  it  is  the  ratio  of  K/Ks^ei  that  determined  whether 
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TABLE  9 


PSOBE  TUBE  EXPERP1ENTAL  RESULTS 
12"  Thick  Fiberqlass  Semple) 


Freq 

(Hz) 

k 

( Rad/n ) 

a 

(Np/n) 

cph 

(n/s) 

125 

14.8 

13.6 

53.1 

250 

20.9 

16.3 

75.2 

500 

29.7 

22.7 

105.8 

1000 

43.3 

21 .7 

145.1 

A 


>-> 


a  material  is  treated  as  a  soft  blanket  or  a  rigid/dence  blanket.  Based  on  the 
appearance,  texture,  and  high  flow  resistivity  of  the  fiberglass  it  was  assumed 
that  the  material  should  be  treated  as  a  rigid/dense  blanket.  [Ref.  4,  pp. 
839-840] 

Once  the  properties  of  the  material  were  characterized,  Eq.  (2.89)  could  be 
used  to  directly  calculate  or  Figure  2.6  and  Figure  2.7  used  to  graphically 
determine  the  magnitude  and  phase  of  the  complex  propagation  constant  b  =  jk  as 
a  function  of  frequency,  from  which  the  complex  wave  number  is  determined. 
Once  these  were  determined  the  characteristic  impedance  could  be  calculated  by 
use  of  Eq.  (2.90).  The  values  determined  from  the  graphs  for  k  and  the 
calculated  values  of  Zc^ar  assuming  isothermal  conditions  in  the  pore  volume 
are  listed  in  Table  10.  The  same  quantities,  assuming  adiabatic  conditions, 
are  listed  in  Table  11. 

As  expected,  the  impedance  magnitude  decreases  with  increasing  frequency. 
The  impedance  phase  calculated  for  both  isothermal  and  adiabatic  conditions  is 
the  same,  a  consequence  of  the  bulk  modulus  being  a  real  quantity  in  the  two 
limiting  cases.  If  isothermal  conditions  were  to  exist,  then  the  impedance 
phase  would  approach  an  angle  of  -45  degrees  (pure  diffusion  of  sound 
pressure) .  Since  the  calculated  phase  values  remain  nearer  to  the  diffusive 
limit,  the  assumption  that  isothermal  conditions  exist  in  the  pore  volume  would 
be  most  nearly  correct. 

Over  the  frequency  range  considered,  the  isothermal  values  for  k  and  a  are 
greater  than  the  adiabatic  values,  which  is  consistent  with  the  fact  that  the 
isothermal  sound  speed  is  less  than  the  adiabatic  sound  speed. 

E.  THE  NEW  TECHNIQUE,  EXPERIMENT  AND  RESULTS 

The  plane  wave  resonator  chosen  to  conduct  all  experiments  was  constructed 
out  of  plexiglass  with  interior  dimensions  as  follows: 


70 


Lenqth  (x-direction)  =  25.4  t  0.1  cn 
Width  (y-direction)  =  19.8  t  0.1  cn 
Heiqht  (z-direction)  =  7.3  ±  0.1  cm. 

This  particular  cavity  was  constructed  for  a  laboratory  experiment  on  standinq 
waves  in  a  rectanqular  room.  The  minimum  wall  thickness  was  1.27  cm  (1/2"). 
The  cavity  was  mounted  on  three  leqs  with  the  source  hole  located  in  the  floor 
at  one  corner  to  allow  the  excitation  of  all  modes.  The  top  cover  of  the  cavity 
is  removable  to  allow  insertion  of  fiberqlass  samples.  When  installed,  the  ton 
cover  seals  the  cavity  and  slides  to  allow  the  receiver  position  to  be  varied. 
A  small  microphone  receptacle  mounted  on  a  slidebar  installed  in  the 

removable  cover  allows  the  pressure  to  be  measured  throughout  the  cavity  with  a 
probe  tube.  All  slidinq  surfaces  were  constructed  with  very  close  tolerances 
and  coated  with  high  vacuum  grease  to  reduce  leakage  and  to  provide  a  tight 
fit  when  the  cavity  cover  was  installed. 

The  receiver  slide  bar  assembly  consisted  of  three  portions:  the  slide  bar, 
the  microphone  receptacle,  and  the  receiver  microphone.  The  slide  bar  fit  into 
a  groove  machined  into  the  interior  surface  of  the  top  cover  and  fit  flush  with 
this  surface.  A  threaded  hole  was  cut  into  the  slidebar  to  allow  mounting  the 
microphone  receptacle.  The  microphone  receptacle  consisted  of  a  hollow 

aluminum  cylinder  with  a  thumb  screw  to  set  and  hold  the  height  of  a  brass 
probe  tube  assembly.  The  probe  tube  assembly  allowed  insertion  of  a  General 
Radio  precision  microphone  flush  with  the  smaller  bore  probe  tube  extension 
which  could  be  inserted  into  the  cavity  interior  without  disrupting  the  sound 
field. 

The  bottom  of  the  cavity  was  narked  with  qrid  lines,  21  in  the  lenqth 

direction,  15  in  the  width  direction,  spaced  1/2"  apart,  producing  336 

individual  grid  points.  The  in-phase  and  quadrature  components  of  pressure 
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TABLE  10 


MATERIAL  PROPERTIES  METHOD  RESULTS 
(Isothermal  Conditions) 

(Rqq  =  15,000  MRS  Ray 1/m) 


Freq. 

(HZ) 

|Zci harl/R0c 

zchar 

(degrees) 

k 

(Rad/m) 

a 

(No/m) 

cph 

(Vs) 

100 

3.66 

-43.5 

6.83 

6.48 

92.0 

125 

3.17 

-43.0 

7.47 

6.96 

105.1 

200 

2.64 

-42.0 

10.08 

9.08 

124.7 

250 

2.32 

-41.5 

11.17 

9.89 

140.6 

300 

2.12 

-41.0 

12.33 

10.72 

152.9 

400 

1.83 

-39.0 

14.66 

11.82 

171.4 

500 

1.66 

-37.5 

16.94 

13.00 

185.5 

600 

1.49 

-37.0 

18.37 

13.84 

205.2 

700 

1.40 

-35.0 

21.26 

14.88 

206.9 

800 

1.32 

-33.5 

22.63 

14.98 

222.1 

900 

1.24 

-32.0 

24.46 

15.28 

231.2 

1000 

1.20 

-31.0 

26.38 

15.85 

238.2 
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TABLE  11 


MATERIAL  PROPERTIES  METHOD  RESULTS 
(Adiabatic  Conditions) 

(RqC  =  15,000  MKS  Ray 1/m) 


Freq.  ~  char  k 

(Hz)  ' zchar ' '  p0c  ( degrees )  (Rad/m) 


a 

(Np/rn) 


were  measured  at  each  qrid  point  for  the  node  and  sample  selected.  Fiqure  3. II 
shows  the  location  of  the  source  and  the  qrid  numbering  scheme  that  was  applied 
to  the  interior  floor  of  the  cavity. 

Ten  fiberglass  samples  were  cut  from  a  fiberglass  sheet,  with  thicknesses 
ranging  from  0.48  cm  (3/16  inch)  to  4.76  cm  (1  7/8  inch)  the  maximum  thickness 
being  equal  to  the  full  thickness  of  the  sheet.  The  samples  were  cut  in  the 
form  of  rectangular  strips  so  as  to  fit  snuqly  alonq  the  side  wall  opposite  the 
source.  Fiqure  3.11  shows  the  location.  The  strips  extended  from  grid  point 
(0,15)  to  grid  point  (20,15).  The  samples  were  cut  such  that  the  surface 
normal  to  the  flow  velocity  in  the  oriqinal  flow  resistance  experiment  would 
also  be  normal  to  the  particle  velocity  expected  when  the  (010)  mode  was 
excited  (i.e.,  the  "width"  direction  in  Figure  3.11)  and  parallel  to  the 
particle  velocity  when  the  (100)  rode  was  excited  (the  "lenqth"  direction). 
Table  12  lists  the  samples  used  and  their  thicknesses. 

Three  experiments  were  conducted: 

1 .  The  mode  structure  in  the  cavity  was  measured  over  the  frequency  ranqe 
of  interest  for  each  sample.  The  effect  of  introducing  the  fiberglass 
samples  on  the  strength  and  shape  of  the  resonances  was  noted.  Of 
particular  interest  were  the  (100)  and  (010)  modes,  which  were  expected 
to  ne  relatively  distinct. 

2.  The  DC  flow  resistance  of  the  air-filled  fiberqlass  was  then  extracted 
from  the  above  measurements  using  the  dispersion  relation  developed 
in  Chapter  2  for  the  allowed  modes  (Eq.  (2.51)),  followinq  which  the 
connlex  wave  number  and  characteristic  impedance  were  calculated. 

3.  The  pressure  field  of  the  1100)  and  (010)  modes  was  measured  for 
varying  sample  thicknesses  and  compared  to  that  predicted  by  theorv. 

The  equipment  setup  used  to  conduct  all  three  of  these  experiments  is  shown 
in  F  ;  jure  3.12. 
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The  first  cf  the  three  experiments  was  conducted  to  determine  the  modes 
present  in  the  empty  cavity  and  the  effects  of  increasing  sample  thickness  on 
the  shape  and  strength  of  each  resonant  mode.  From  mode  theory,  seven  modes 
were  expected  to  exist  in  the  empty  cavity  for  frequencies  less  than  2000  Hz. 
Of  these  seven  modes,  the  two  lcwest  were  the  [100]  and  [010]  modes,  both  of 
which  occur  below  1000  Hz.  All  seven  modes  were  observed,  with  the  two  lowest 
modes  being  strong  and  relatively  distinct  compared  to  the  others.  This  is 
shown  in  Figure  3.13.  Table  13  lists  the  theoretically  calculated  and 
experimentally  measured  peak  frequencies  of  the  modes  present  in  the  cavity  up 


There  was  some  concern  that  as  thicker  samples  were  introduced  the 
resonance  peaks  would  overlap,  makinq  the  measurement  of  the  half-power  widths 
unreliable.  This  condition  set  the  upper  limit  on  usable  sample  thickness  to  1 
7/8”  for  the  [100]  mode  (sample  #10)  and  11/16"  for  the  [010]  mode  (sample 
#5). 

For  all  observed  modes,  the  center  frequency  down-shifted  as  the  sample 
thickness  was  increased,  while  the  0  decreased.  The  [100]  mode  has  an  empty 
cavity  peak  frequency  and  Q  of  680.2  Hz  and  147.9,  respectively.  Inserting  the 
1  7/8"  (sample  #10)  fiberglass  strip  down-shifted  its  center  frequency  and 
lowered  the  0  to  570.0  Hz  and  3.2,  respectively,  for  a  total  frequency  change 
of  110.2  Hz.  The  [010]  mode  has  a  center  frequency  of  870.2  Hz  and  a  0  of  57.6 
in  the  empty  cavity,  but  when  the  11/16”  (sample  #5)  fiberglass  strip  was 
inserted  the  peak  frequency  only  down-shifted  17  Hz  to  853.2  Hz  with  a  Q  of 
21.1.  For  both  mode  measurements  the  receiver  probe  tube  remained  flush  with 
the  ceiling  of  the  cavity  at  grid  point  (20,0).  See  Fiqures  3.14  and  3.15 
Several  observations  from  this  experiment  were  made: 

(1)  The  0  and  the  center  frequency  of  both  modes  decreased  with  increasing 
sample  thickness. 

(2)  Compared  to  the  [010]  mode,  the  [100]  mode,  with  particle  velocity 
parallel  to  the  fiberglass  surface,  allowed  thicker  samples  to  be 
introduced  before  the  resonance  curve  spread  such  that  the  Q  could  no 
longer  be  reliably  determined. 

(3)  For  all  sample  thicknesses,  the  0  of  the  resonance  curves  was 
larger  with  the  receiver  positioned  above  the  source  (grid  point  (0,0)) 
than  in  the  opposite  corner  (grid  point  (20,0)).  This  was  assumed 
to  be  due  to  near  field  effects  of  the  source  on  the  receiver. 
Consequently,  only  the  resonance  data  taken  with  the  receiver  at  grid 
point  (20,0)  is  reported  here. 
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Sample  # 
1 
2 

3 

4 

5 

6 

7 

8 

9 

10 


TABLE  12 

FIBERGLASS  SAMPLE  THICKNESS 


! 

1 


Thickness 
( inches) 

Thickness 

(cm) 

3/16  ±  1/32 

0.48  ±  0.01 

1/4 

0.64 

3/8 

0.95 

1/2 

1.27 

11/16 

1.75 

7/8 

2.22 

1 

2.54 

1  1/4 

3.18 

1  1/2 

3.81 

1  7/8 

4.76 
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Figure  3.12  Experimental  Setup  for  the  New  Method 
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TABLE  13 


RESONANCE  MODES  OF  THE  EMPTY  CAVITY 
(Theoretical  and  Measured) 


Mode  # 

Resonant 

Frequency 

(calculated) 

Resonant 

Frequency 

(measured) 

100 

679 

£ 

4 

676  £  1 

010 

871 

♦ 

5 

868 

110 

1105 

9 

1099 

200 

1358 

♦ 

12 

1354 

210 

1614 

14 

1608 

020 

1742 

♦ 

19 

1738 

120 

1870 

♦ 

23 

1865 

FW.QO***.'* 


Fiqure  3.13  Mode  Structure  of  the  Qnpty  Cavity 


'Q 


(4)  Overlap  of  the  [110]  mode  with  the  [010]  node  was  evident  for 
fiberqlass  thickness  of  1/2"  (sample  #4)  and  greater.  Overlap  of 

the  [100]  node  with  the  [010]  node  was  not  significant  until  a 
thickness  of  1  1/4"  (sample  #8)  of  fiberglass  was  introduced.  See 
Fiqures  3.14  and  3.15. 

2 .  *teasured  fc  and  Q  of  the  [100]  and  [010]  Ntodes  for  Various  Fiberqlass 
Thicknesses 

The  second  experiment  consisted  of  measuring  the  center  frequency  and  0 
of  the  [100]  and  [010]  modes  as  a  function  of  sample  thickness.  The  receiver 
was  placed  in  the  ceiling  of  the  cavity  at  grid  point  (20,0)  at  the  opposite 
end  of  the  long  wall  from  the  source.  This  location  was  chosen  to  minimize 
direct  acoustic  cross-talk  between  the  source  and  receiver.  For  the  [100] 
mode,  each  of  the  10  samples  was  inserted  and  the  center  frequency  and  Q  were 
measured  from  the  resonance  curve.  This  sample  procedure  was  repeated 
for  the  [010]  mode  except  only  samples  1  through  5  were  used.  The  center 
frequency  and  half-pcwer  frequencies  were  mesured  to  within  a  tenth  of  a  Hertz 
by  readinq  the  frequency  of  the  1  to  1.5  fUz  Local  Oscillator  output  of  the 
'P3580A  Spectrum  Analyzer  on  a  HP5300B  Frequency  Counter.  Settinq  the  scale  to 
L  IB/div  on  the  Spectrum  Analyzer  allowed  the  center  frequency  and  half-power 
full  width  of  the  resonance  curves  to  be  measured.  Tables  14  and  15  list  the 
experimentally  determined  values  of  fc  and  0  for  the  [100]  and  [010]  modes. 

Following  the  experimental  measurements,  theoretical  values  of  fc  and 
»•  for  the  [100]  and  [010]  modes  as  a  function  of  sample  thickness  were 
ralculated  frcm  the  complex  roots  of  the  dispersion  relation  (Eq.  (2.51))  for 
various  values  of  DC  flew  resistance  using  a  numerical  analysis  technique  known 
as  .Sard’s  method  [Ref.  19,  pp.  169-170].  A  brief  description  of  this  method  is 
liven  m  Appendix  A.  However,  before  a  comparison  between  theory  and 
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Figure  3.14  Mode  Structure  of  Cavity  With  a  11/16"  Fiberglass  Strip 
(sample  #5) 
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Figure  3.15  Mode  Structure  of  the  Cavity  With  1  7/8"  Fiberglass  Strip 
(sample  #10) 


81 


«xner inert  )  '*?  marie ,  and  the  nest  fittinn  value  >f  'jr  t  i  -w  rps;;Mr  « 

Jetemined,  :*•  was  oe<essar/  t-  j^rrect  me  ■*•  the  irner  tor  *wnr>t  /  .  =»■;/  / 
l<>sses.  vjur  :*-s  >f  *rrf  v  :av  1  r  /  ; usses  are  .■ons  i)ere<l  i"  •■•'.*»  oe*' 

suhsect  i  m  . 

j .  v»,r  >*s  t  '.jjh s  i'  - *ie  tmry  ~av  1 t  / 

I-1  lerivinq  *:ne  lisoersi  m  relation  for  trie  new  r*»sr/nan<;*»  ’•*•  ••  .  ;.»*  .• 

was  assumed  ••hat  *■  re  r\ iy  losses  wh  1  oh  jccur  ir  tr*  oiane  wave  r»>s  .ma *■  ,r  ao- 
lue  to  v  1  scouS  anti  thenal  1'tsses  wittun  the  fijid  sat  ..rated  -»ir'»,s  a.u  t . 
Although  these  are  'he  largest  ^mr.tts  of  loss  in  the  cavity,  other  1  »s 
mechanisms  must  me  accounted  tor.  To  account  for  these,  several  assu*n>i  i- ms 
were  made : 

II)  trie  losses  in  toe  empty  cavity  remain  the  sane  ever,  trviuqh  a 
nort  ion  of  the  volijne  and  wall  surface  is  lost  lue  to  the  i  ntrojuct  i  on 
of  a  fibemlass  sarnie  'less  than  10%), 
i  2 1  the  shift  of  resonance  f  re*;uency  lue  to  the  introduction  if  i  n.ner- 
3 lass  sample  does  not  effect  the  j  of  the  «spty  cavity  and, 

'3)  each  loss  mechanism  ;an  he  treatetl  as  act  1  rrj  l mlenendent ly . 

,ene  rally,  stater),  it  was  assumed  that , 

Total  lasses  m  the  »  losses  m  the  ♦  1  »ses  Jue  to 

fiberglass  filler)  empty  :avity  only  to  fiber- 

cavity  class 

which  can  he  written  in  terns  of  juality  factors  -as 

4 -  *  i-  ♦  ,  (3.6) 

OpTT  ^F 

where  l>jr7r  13  0  of  the  cavity  lue  to  all  loes  mechanisms,  Op-  is  the  emoty 

cavity  0  and  Op  is  the  Q  due  to  viscous  and  thermal  losses  within  the 
fiberglass.  Eq.  (3.6)  was  used  to  crrpare  the  measured  values  of  Oj^  to  the 
calculated  values  of  Op  for  the  (1001  and  (010)  modes  by  setting  the  values  of 
Op  equal  to  those  experimentally  determined  for  the  empty  cavity. 
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TABLE  14 


1100) 

'COE  '.'ENTER 

FREQUENCY  AND  0 

VALUES 

Thickness 

fc 

U 

fu 

Sample  • 

(cm) 

(Hz) 

(HZ) 

(HZ) 

Onptv 

0.0 

680.2  ♦.  0. 

1  678.2  t  0.1 

682.8  t  0.1 

14^.9 

1 

0.48 

669.5 

663.9 

675.3 

58.7 

2 

0.64 

667.5 

661.3 

674.1 

52.2 

1 

0.95 

658.6 

648.8 

668.1 

34.1 

4 

1.27 

652.1 

640.1 

664.2 

27.1 

5 

1.75 

646.3 

628.9 

659.9 

20.9 

6 

2.22 

623.6 

587.0 

650.5 

9.8 

7 

2.54 

614.1 

567.4 

647.6 

7.7 

8 

3.18 

592.3 

531.9 

629.7 

6.1 

9 

3.81 

571.6 

506.6 

624.3 

4.9 

10 

4.76 

570.0 

495.4 

672.6 

3.2 

83 


j^UL 


TABLE  15 

[010]  fCDE  CENTER  FREQUENCY  AND  0  VALUES 


Thickness 

Ec 

E* 

£u 

Sample  * 

(cm) 

(Hz) 

(Hz) 

u 

(Hz) 

Q 

Empty 

0.0 

870.2  t  0.1 

863.0  i  0.1 

878.1  t  0.1 

57.6 

1 

0.48 

865.1 

854.5 

874.4 

43.5 

2 

0.64 

864.4 

852.0 

876.2 

35.7 

3 

0.95 

858.8 

846.4 

872.3 

33.2 

4 

1.27 

856.5 

841.1 

873.1 

26.8 

5 

1.75 

853.2 

834.4 

874.9 

21.1 
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An  attempt  to  identify  the  major  source  of  loss  in  the  empty  cavity  was 
made,  as  tnis  in  turn  determined  the  value  of  Og  to  be  used  in  correcting  for 
losses  other  than  in  the  fiberglass.  The  loss  mechanisms  in  the  empty  cavity 
were  assumed  to  be  due  to  3  mechanisms: 

—  transmission  of  energy  through  the  plexiglass  walls  of  the  cavity, 

—  viscous  ard  thermal  losses  at  the  interior  surfaces  of  the  cavity, 

—  leakage  of  energy  due  to  incomplete  sealing  of  the  top  cover  and  sliding 
section  of  the  cavity  and  slide  bar  assembly. 

The  transmission  of  energy  through  the  walls  was  of  concern  mainly  since  the 

cavity  was  not  designed  specifically  for  this  experiment.  For  normal 

incidence,  the  transmission  loss  through  the  walls  of  the  cavity  was  calculated 

to  be  32  dB  for  the  minimum  wall  thickness  (1.27  cm),  [Ref.  17,  p.  128].  Based 

on  the  maximum  wall  thickness  (2.54  cm)  the  lowest  expected  coincidence 

freguency  for  a  free  plexiglass  plate  was  calculated  to  be  1282  Hz.  For  a 

plate  fixed  about  its  edges  this  freguency  is  higher.  Therefore,  transmission 

by  coupling  to  flexural  waves  in  the  walls  can  be  ignored.  The  following 


eguation  was  used  to  estimate  the  0  resulting  from  each  source  of  loss: 

_  2  n  Estored 


Elost/cycle 


(3.7) 


where 


Estored  *  cV 


(3.8) 


and  e  is  the  average  energy  density,  and  V  is  the  interior  volume  of  the 
cavity.  The  energy  lost  per  cycle  can  be  written 


1 

Elost/cycle  2  ~4  e  TI  f  ' 


(3.9) 


where  A  is  the  surface  area,  c  is  the  speed  of  sound  in  the  fluid,  Tj  the 
power  transmission  coefficient  and  f  is  the  freguency.  Inserting  Eg.  (3.8) 
and  (3.9)  into  Eg.  (3.7)  results  in 


.WA.Vv'' 


*» 


Qtrans 


<  3.10  l 


9nV 

~  XTjA  * 

For  loss  due  to  transnission  throuqht  the  cavity  walls,  assuming  1000 
Hz  as  the  worst  case  (X  *  0.343  m) ,  the  value  of  Otrans  was  calculated  to  ne 
2700.  It  was  concluded  that  transmission  throuqh  the  walls  was  not 
significant. 

Next  considered  were  the  effects  of  the  viscous  and  thermal  losses 
occurrinq  at  the  cavity  surfaces.  Aqain  Eq.  (3.7)  was  used  to  determine  3surf» 
with  Egtored  ceroaifuncj  the  same  as  in  Eq.  (3.8).  The  energy  lost  per 
cycle  due  to  viscous  effects  can  be  written  approximately  as 

E lost/cycle  2  ^visc  1.11' 

where  6visc  the  viscous  penetration  depth  in  air.  The  energy  lost  per 

cycle  due  to  thermal  effect  can  be  written  approximately  as 

Elost/cycle  2  **  ( T~1 ) ^th  (3.12) 

where  y  is  the  ratio  of  specific  heats  Cp/cv  (*  1.4)  for  the  air  anti  is 
the  thermal  penetration  depths  in  air.  Since  the  viscous  and  thermal 
penetration  depths  in  air  are  approximately  equal,  an  expression  for  the 
energy  lost  per  cycle  due  to  both  viscous  and  thermal  effects  is 

E lost/cycle  2  ^isc  y  '  '3.13) 

In  other  words, 

«visc  +  ( r-i 3  2  *«visc  ♦  (3-14' 

Substituting  Eq.  (3.8)  and  (3.13)  into  Eq.  (3.7)  results  in 


Qsurf 


2irV 

Advise  Y 


( 3. IS) 


Choosing  a  frequency  of  100  Hz  for  the  worst  case  (greatest  viscous  and  thermal 


penetration  depths),  with  the  air  temperature  in  the  cavity  assumed  to  be  20#C, 
the  resulting  0gurf  was  450.  Therefore,  viscous  and  thermal  wall  losses  are 
negligible. 


Lastly,  the  loss  flue  to  leakage  oast  cavity  slideinq  surfaces  was 
considered.  Since  this  loss  could  not  he  directly  calculated,  a  measure  >t 
its  significance  was  estimated  Py  pressing  down  on  the  too  cover  cf  the  cavity 
at  the  wall  edges  while  retakinq  the  [100!  mode  resonance  curve  lata  for  each 
sarnie.  The  empty  cavity  0  increased  from  147  to  262,  with  a  1 . -)  Hz  inshift  m 
the  footer  freguency.  From  this  experiment,  a  value  if  j  for  leaxa>je  was 
estimated  to  he  approximately  280.  This  indicated  that  leakage  was  the  naor 
source  of  loss  in  the  empty  cavity  ty  at  least  a  factor  of  two.  rh*»  leak.t  je 
problem  could  not  he  corrected  without  making  maior  changes  in  the  caviw 
design.  No  attempt  was  made  to  account  for  the  upshift  in  the  center  frwjuencv 
that  occurred  wtule  pressing  down  x\  the  cavity  top  (the  frequency  ipshift  was 
xi  the  order  of  1  H2  for  each  sample,  x4ule  the  frequency  downshift  lue  *  > 

introducing  the  fiherglass  was  on  the  average  of  10  to  12  Hz  drop  ner  sample). 

4.  Determinat  ion  of  Rrjr  from  fc  and  J 

For  the  (1001  mode,  the  iterative  process  used  to  extract  the  '*  f  1  w 
resistance  of  the  fiherglass  from  the  dispersion  relation  ( F.g.  '  2 .  S  l  >  >  was 
started  hy  choosing  equal  to  V»,200  *ecS  Ray  1  'm  (the  mean  flow  resistance 
letermined  in  the  large  sample  hol<Jer  with  particle  velocity  mrmal  to  the 
x  \  r-f  iher'jlass  surface).  After  entering  the  fluid  ami  cavity  parameters  (see 
Appendix  A  for  ;>arameters  ami  values  used),  the  roots  of  the  dispersion 
relation  for  *  were  then  determined  for  each  sample  usimj  the  measure)  values 
H  f,.  xml  J  m  an  initial  starting  point.  The  resulting  nnts  for  each  sample 
were  numerically  computed  until  the  sum  of  the  magnitudes  of  the  real  ami 
imaqinery  :»art*  of  the  dispersion  relatrxi  (Fg.  2.61)  was  less  than  1  x  lf)“b. 
Hy  this  means,  -alculared  values  of  vs  fc  ware  obtained  for  a  particular 
value  of  Rfy. 

The  calculated  ami  measured  values  of  0  vs  tv.  were  compared  in  two 
ways.  Figure  1.16  shows  vs  fr.  Here  the  calculated  values  of  0^  have 
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been  corrected  for  emoty  cavity  losses  using  Eq.  (3.6).  Fiqure  3.17  shows  iy 
vs  fc.  Here  the  measured  values  of  have  been  corrected  for  empty  cavity 

losses  using 

j_  _  1 _  _  1 

Op  $rcrr  ^e* 

If  the  calculated  curve  plotted  above  the  measured  curve,  this 
indicated  too  hiqh  a  value  of  Rqc  and  a  lcwer  value  of  was  chosen.  If  the 
curve  plotted  below  the  measured  curve,  this  indicated  too  lew  a  value  of  R^ 
and  a  larger  value  of  Rq^  was  chosen.  This  process  was  continued  until  the 
measured  curve  was  bounded  between  two  values  of  Rq£ •  Once  bounded,  smaller 
increments  were  chosen  for  Rpc  and  bh«  process  was  continued  until  a  Pest  fit 
(graohical ly)  was  obtained  between  the  calculated  and  measured  curves 
of  Op  or  Ojot  vs  fc.  For  the  (100)  mode,  the  best  fit  value  for  R^  was 

15,000  rus  Rayl/m,  which  was  four  times  less  than  that  measured  experimentally. 

Two  possible  causes  were  investigated  to  explain  this  discrepency.  The 
first  possible  cause  is  the  10%  compression  of  the  samples  observed  in  the 
original  PC  flow  resistance  experiment  (large  sample  holder).  \ssumnq  the 
total  pressure  drop  across  the  test  sarnie  was  independent  of  the  compressed 
thickness,  the  extracted  value  of  DC  flew  resistance  per  unit  length  should  be 
corrected  trm  59,200  fflCS  Rayl/m  to  53,000  l*S  Rayl/m.  The  second  possible 

cause  is  that  the  flow  resistance  of  the  fiberglas  is  not  isotrootc,  and  that 

it  is  greater  for  flow  perpendicular  to  the  sheet  faces.  This  might  be 
expected,  as  there  is  visual  evidence  that  the  fiberglass  sheets  have  been 
compressed.  In  the  first  flow  resistance  measurement,  the  fiberglass  sample 
was  oriented  so  that  the  fiberqlass  sheet  surface  was  normal  to  the  flow 
velocity.  However,  for  the  (1001  mode,  the  particle  velocity  is  parallel  to 
the  fiberglass  sheet  surface.  Flow  resistance  tests  were  again  conducted, 
using  a  much  analler  sample  holder,  with  samples  cut  so  that  the  anisotropy  in 
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*a  Center  Frequency 


CENTER  FREQUENCY 


vs  Center  Frequency 


Rq£  could  be  measured.  As  indicated  earlier  in  this  chapter,  the  flow 
resistance  normal  to  the  sheet  surface  now  measured  approximately  three  times 
greater  than  that  parallel  to  the  sheet  surface.  The  value  measured  with  flow 
parallel  to  the  surface  was  15,200  fSCS  Rayl/m.  With  this  value  of  Rqq,  the 
roots  of  the  dispersion  relation  were  recalculated,  Opox  and  calculated  and 
plotted.  These  values  fell  within  8%  of  measured  values  for  Oror  vs  and 
within  12%  for  Op  vs  fc  for  all  samples. 

The  method  of  extracting  from  measured  values  of  fc  and  0  was  then 
applied  to  the  (010)  mode,  which  has  particle  velocity  normal  to  the  sheet 
surface.  Startinq  with  a  value  of  34,000  MKS  Rayl/m  for  Rpc  ( the  value 
measured  in  the  small  sample  holder),  values  of  0  vs  fc  were  calculated  for 
various  values  of  Rq-  and  were  compared  with  measured  values  of  Q  vs  fc.  The 
measured  values  of  0  vs  fc  are  founded  by  the  calculated  curves  for 
RqC  *  5,000  MCS  Rayl/m  and  »  20,000  MKS  Rayl/m,  and  are  consistent  with  a 
value  of  15,000  rt(S  Rayl/m  for  the  same  as  for  the  [100)  mode.  This  is 
shown  in  Fiqures  3.18  and  3.19.  The  reason  for  the  discrepency  between  these 
boundinq  values  and  the  measured  value  of  Rq^  for  flow  normal  to  the  fiberglass 
sheet  face  is  not  known.  It  may  be  due  to  mode  overlap  of  both  the  [110)  and 
[100)  modes  on  to  the  (010)  mode.  Fiqures  3.20  and  3.21  show  plots  of  the 
resonance  center  frequency  for  the  [100)  and  [0101  modes,  respectively,  as  a 
function  of  sample  thickness  for  the  same  values  of  Rqq  as  in  Fiqures  3.18 
and  3.19. 

Having  determined  that  the  best  fit  value  of  the  EX?  flow  resistance  in 
the  porous  material  was  approximately  15,000  MKS  Rayl/m,  the  complex  wave 
number  and  characteristic  impedance  were  calcualted  usinq  Eq.  (2.56)  and  Eq. 
(2.57)  and  the  roots  of  the  dispersion  relation,  (Eq.  (2.51)).  Table  16  lists 
the  results  for  the  (100)  mode,  while  Table  17  lists  those  for  the  [010)  mode. 
The  data  for  both  modes  is  graphically  presented  in  Fiqures  3.22  through  3.26. 
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Sensitivity  of  Roots  of  the  Dispersion  Relation  to  Variation  n 
Parameters 

In  order  to  determine  the  sensitivity  of  the  roots  of  the  dispersion 
relation  to  variations  in  fluid  and  porous  material  physical  properties,  the 
roots  were  aqain  calculated  for  the  [100]  mode  with  each  of  the  input 
parameters  chanqed  toy  a  small  amount,  one  at  a  time,  and  the  character ist 10 
impedance  and  complex  wave  number  were  recalculated. 

It  was  found  that  there  was  less  than  a  1%  chanqe  m  •<  and  z  .„ar 
when  porosity  was  decreased  to  0.95  (5%  chanqe),  when  bulk  fluid  sound  speed 
was  chanqed  to  343  iVs  (0.6%  chanqe  or  2  m/s),  and  when  bulk  fluid  lensitv  was 
changed  to  1.23  Kq/m3  (1.6%  chanqe  or  0.02  Kq/m3).  However,  an  11%  change  m 
and  Zcpar  was  found  when  the  DC  flow  resistance  was  varied  oy  ♦  s.ooo  MK> 
Rayl/m  about  15,000  MCS  Rayl/m.  Because  the  uncertainty  in  the  measure-1  lata 
is  so  much  greater,  the  sensitivity  of  the  best  fit  value  of  R^,  and  nence  «. 
and  Zchar#  to  slight  uncertainties  in  other  physical  properties  is  negi :  nr 'Le. 


6 .  Modal  Pressure  Field  in  a  Rectangular  Cavity  with  Porous  Vroust  i 
Material  . 

To  theoretically  determine  the  (100)  and  (010)  modal  pressure 
distribution  in  the  cavity  for  each  fibecqlass  sample  requires  the  values  ->t 
k^z  and  k2z  be  determined  from  the  roots  of  the  dispersion  relation  for 
for  the  best  fit  value  of  Rq£.  Then,  usinq  Eq.  (2.63)  for  the  bulk  fluid  and 
Eq.  (2.64)  for  the  porous  material,  the  normalized  r’aqmtude  and  phase  the 
pressure  at  each  qrid  point  can  he  determined.  See  Appendix  B  for  a  more 

detailed  discussion.  For  the  (100)  mode,  Figure  3.27  shows  graphically  the 


results  for  sample  #1  ( 3/16" )  and  Figure  3.28  shows  the  results  for  sample  110 


(1  7/8").  For  the  (010]  mode,  Fiqure  3.31  shows  graphically  the  results  for 
sample  #1  (3/16")  and  Figure  3.32  shows  the  results  for  sample  #5  (11  lb"). 
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Figure  3.21  1010]  Mode  fc  vs  Sample  Thickness 
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Figure  3.24  k  vs  f  For  New  Method  ( [100}  and  [010]  Modes) 
Absorption  Coeff.  vs  Frequency 

(NEW  TECHNIQUE— Cl 003  l  C0103) 
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Phase  Speed  vs  Frequency 

(NEW  TECHNIQUE— n 003  &  C0103) 


Figure  3.26  Cp^  vs  £  For  New  Method  ([100)  and  [010]  Modes) 
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TABLE  16 


ZChar 

(Roc 

k  RESULTS 
■  15,200 

FOR  (100)  NODE 
f«S  Ray 1/m) 

Sample  1 

tc 

(HZ) 

k 

(Rad/m) 

a 

(Np/ta) 

l^charl/ Qoc 

$ 2  u 
£char 

(deqrees) 

°ph 

(m/s) 

1 

669.4 

±0.1  21.54 

14.53 

1.48 

-37.2 

195.2 

2 

665.4 

21.47 

14.47 

1.48 

-37.2 

194.8 

3 

657.7 

21.32 

14.35 

1.49 

-37.1 

193.8 

4 

649.8 

21.17 

14.21 

1.50 

-37.1 

192.9 

5 

637.6 

20.95 

13.98 

1.52 

-37.0 

191.2 

6 

625.3 

20.74 

13.71 

1.53 

-36.8 

189.4 

7 

617.1 

20.61 

13.52 

1.54 

-36.6 

188.1 

8 

601.5 

20.38 

13.10 

1.56 

-36.2 

185.4 

9 

587.4 

20.19 

12.66 

1.58 

-35.6 

182.8 

10 

569.5 

20.03 

11.96 

1.61 

-34.4 

178.6 
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TABLE  17 


Zchar  AND  k 

'  RESULTS 

TOR  [010]  MODE 

( &DC  * 

15,000 

MKS  Rayl/tn) 

fc 

(Hz) 

k 

(Rad/m) 

& 

(Np/m) 

1 2*char  1  /  Poc 

£char 

(deqrees) 

cph 

(n/s) 

865.0  i  0.1 

25.76 

15.67 

1.30 

-34.8 

2 13.5 

862.9 

25.72 

15.63 

1.30 

-34.8 

210.8 

858.5 

25.64 

15.55 

1.31 

-34.8 

210.4 

854.3 

25.57 

15.46 

1.31 

-34.8 

209.9 

848.3 

25.47 

15.29 

1.32 

-34.7 

209.3 

101 


These  plots  present  the  node's  nodal  pressure  field  in  phasor  notation.  The 
naqnitude  is  represented  by  the  size  of  the  arrow  and  the  ohase  anqle  by  the 
angular  rotation  of  the  arrow  in  degrees  from  the  horizontal  (i.e.,  zero 
degrees  of  ohase  angle  is  a  phasor  pointing  horizontally  to  the  right). 

Experimentally,  the  pressure  field  was  measured  by  inserting  each 
sample  into  the  cavity  and  exciting  the  desired  mode  at  its  center  frequency 
for  that  sample.  By  programing  a  HP-300  series  oomputer  to  take  data  from  die 
SR  530  Dual  Channel  Lock-In  Amplifier  (in  the  X,  Y  mode),  the  in-phase  and 
quadrature  components  of  pressure  were  measured  at  each  grid  point.  The 
normalized  magnitude  and  phase  of  the  pressure  was  then  plotted.  For  the  [100] 
mode,  Figure  3.29  displays  the  results  for  sample  #1  (3/16")  and  Figure  3.30 
displays  the  results  for  sample  #10  (1  7/8").  The  results  for  the  [010]  mode 
are  displayed  in  Figures  3.33  and  3.34  for  sample  #1  (3/16")  and  sample  #5 
(11/16"),  respectively. 

In  comparing  the  theoretical  and  experimental  pressure  field  plots 
several  observations  can  be  made: 

(1)  The  effects  of  cavity  leakaqe  at  the  slidinq  surfaces  is  evident 
when  the  experimental  plot  (Fiqure  3.29)  is  compared  to  the  theoretical 
plot  (Fiqure  3.27)  of  the  (100)  mode  for  sample  #1  (3/16")  by  noting 
that  the  pressure  phasors  alonq  the  boundaries  not  in  contact  with  the 
fiberglass  strip  in  Figure  3.29  ex i bit  a  small  increase  in  phase  anqle 
(i.e.,  they  are  not  horizontal  as  they  would  be  if  no  losses  existed  at 
the  boundaries)  which  is  not  present  in  Fiqure  3.27.  This  effect  nay 
also  be  seen  for  the  (010)  mode  in  Fiqure  3.31  and  Fiqure  3.33. 

(2)  The  effects  of  the  porous  material  on  the  pressure  distribution  in 
the  hulk  fluid  can  be  seen  to  "qrow"  with  increasing  sample  thickness  by 
observing  the  relative  increase  in  the  phase  anqle  away  from  the 
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[100]  Mode,  3/16"  (  sample  #1  ) 
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air-fiberqlass  boundary  fron  the  theoretical  plot  of  the  3/16"  (sample 
#1)  fiberqlass  strip  compared  to  the  theoretical  olot  of  the  1  7/8" 
(sample  #10)  fiberqlass  strip.  This  can  be  seen  in  Figures  3.27  and 
3.28,  respectively.  This  effect  is  not  seen  in  the  theoretical  plots  of 
the  (010)  mode  for  sample  #1  (Fiqure  3.31)  and  for  sample  #5  (Figure 

3.32) 

(3)  The  effects  of  mode  overlap  become  more  evident  as  sample  thickness  is 
increased  by  comparing  the  measured  and  calculated  modal  pressure  fields 
for  the  same  sample  thickness.  For  example,  compare  Fiqure  3.29  (sample 
#1)  to  Fiqure  3.30  (sample  #10)  for  the  [100]  mode  and  Fiqure  3.33 
(sample  #1)  to  Figure  3.34  (sample  #5)  for  the  [010]  mode. 

(4)  The  pressure  distribution  at  the  air-fiberqlass  boundary  shows  the 
effect  of  the  relative  size  and  sign  of  k\z  and  k2z  by  noting  that  for 
the  thickest  samples  the  pressure  phasor's  anqle  increases  as  the 
interface  is  approached.  For  the  [100]  mode  compare  the  theoretical  plot 
for  sample  #1  (Figure  3.27)  to  the  theoretical  plot  for  sample  #10. 
(Figure  3.28)  and  for  the  [010]  mode  compare  the  theoretical  plot  for 
sample  #1  (Figure  3.31)  to  the  theoretical  plot  for  sample  #5  (Fiqure 

3.32) . 

The  primary  reason  for  comparing  the  calculated  and  measured  pressure 
distributions  was  to  verify  which  node  was  excited  and  to  check  its  purity.  The 
determination  of  whether  the  expected  mode  had  been  excited  was  made  by 
comparing  the  calculated  and  measured  pressure  field  plots  for  the  thinnest 
fiberglass  strip  (3/16"),  For  the  [100]  mode  the  nodal  plane  runs  normal  to  the 
length  direction  from  qrid  point  (10,0)  to  qrid  point  (10,15).  Travelling 
across  the  nodal  plane,  a  180  degree  phase  shift  in  the  pressure  phasor  is 
observed,  with  the  pressure  magnitude  reachinq  a  maximum  at  the  cavity 
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boundaries  normal  to  the  lenqth  direction.  Likewise,  for  the  [0 10 j  mode  the 
nodal  plane  is  normal  to  the  width  direction  fallinq  half-wav  between  the 
boundaries  parallel  to  the  lenqth  direction.  The  pressure  phase  also  chanqes 
180  deqrees  in  crossinq  this  nodal  plane,  with  its  maqnitude  reachinq  a  maximum 
at  the  fiberglass  boundary  and  opposite  wall. 

The  purity  of  these  modes  is  retained  only  with  the  thinner  samples  and 
becomes  contaminated  as  thicker  samples  are  introduced.  This  is  clearly  seen  by 
comparinq  Fiqures  3.29  and  3.30. 
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IV.  COMPARISON  OF  RESULTS  AND  CONCLUSION 


A.  COMPARISON  OF  TOE  COMPLEX  fcftVE  NUMBER 

The  complex  wave  number  may  be  written  in  terns  of  its  real  and  imaqinery 
parts  as  the  wave  number  k  and  the  absorption  coefficient  a  ,  resoectively. 
Figure  4.1  shows  the  wave  number  k  as  a  function  of  frequency  for  each  method. 
Values  of  k  determined  by  the  new  method  fall  between  those  found  by  Yaniv's 
method  (the  upper  bound)  and  by  Beranek's  Material  Properties  method  in  the 
isothermal  limit  (the  lower  bound).  For  example,  at  600  Hz  the  value  of  the  wave 
number  determined  by  the  new  method  was  within  6%  of  the  value  determined  by 
Yaniv's  method.  The  values  of  k  from  the  probe  tube  experiment  are  about  twice 
those  determined  by  the  new  method.  This  is  to  be  expected,  since  the  particle 
velocity  in  the  probe  tube  experiment  was  normal  to  the  air-fiberglass  boundary, 
i.e.,  was  in  the  direction  of  higher  measured  flow  resistance.  (Note  that 
although  the  1010]  mode  also  has  particle  velocity  normal  to  the  fiberglass 
surface,  the  best  fit  value  of  Rqq  (15,000  MCS  Ray 1/m)  was  used  to  calculate  both 
the  complex  wave  number  and  characteristic  impedance). 

The  abosorption  coefficient  as  a  function  of  frequency  is  displayed  in  Fiqure 
4.2  for  each  method.  The  values  determined  by  the  new  method  fall  about  those 
obtained  by  Beranek's  Material  Properties  method  (in  the  isothermal  limit).  For 
example,  at  600  Hz  the  value  of  the  absorption  coefficient  determined  from  the 
new  method  falls  with  17%  of  the  value  determined  by  Yaniv's  method.  Again,  the 
values  determined  usinq  the  probe  tube  method  were  about  twice  those  measured  by 
the  new  method.  The  cause  of  the  dip  in  the  values  of  the  absorption  coefficient 
determined  by  Yaniv's  method  is  not  known. 

Figure  4.3  displays  the  phase  speed  in  the  porous  material  as  a  function  of 
frequency  for  each  method.  Values  determined  by  the  new  method  are  bounded  by 


FREQUENCY 


those  found  by  Yaniv's  method  and  by  Beranek's  Material  Properties  method  (in 
the  isothermal  limit).  For  example,  at  600  Hz  the  phase  speed  determined  by  the 
new  method  is  within  7%  of  the  value  determined  by  Yaniv's  method.  The  values 
of  phase  speed  determined  by  the  probe  tube  method  are  approximately  one-half  of 
those  found  by  the  new  method  for  the  same  reason  discussed  above. 

B.  COMPARISON  OF  TOE  CHARACTERISTIC  IMPEDANCE 

The  normalized  characteristic  impedance  as  a  function  of  frequency  for  all 
methods  is  displayed  in  Figure  4.4.  Values  determined  by  the  new  method  are 
bounded  by  those  from  the  adiabatic  and  the  isothermal  material  properties 
method.  Values  from  the  new  method  plot  almost  directly  on  too  of  the  Biot 
theory  in  the  diffusion  limit.  For  example,  at  600  Hz  the  value  determined  by 
the  new  method  fell  within  26%  of  that  determined  by  Yaniv's  method  and  within 
20%  of  the  impedance  tube  value. 

Figure  4.5  displays  the  impedance  phase  as  a  function  of  frequency  for  each 
method.  The  values  determined  by  the  new  method  fall  about  the  values 
determined  by  the  material  properties  method.  The  low  frequency  values  (<  200 
Hz)  determined  by  the  impedance  tube  method  are  thought  to  be  in  error  due  to 
the  difficulties  in  measuring  the  distance  to  the  first  node  in  that  method.  It 
is  not  known  why  the  values  of  the  phase  extracted  using  the  new  method  frcm  the 
[100]  mode  data  increase  for  decreasing  frequency,  since  in  the  low  frequency 
limit,  the  phase  should  approach  -45  degrees  (pure  diffusion  of  sound  pressure). 

C.  CONCLUSION 

The  ability  to  determine  a  material's  DC  flow  resistance  from  parameters  of 
the  resonance  curve  in  a  plane  wave  resonator  was  verified.  Additionally,  the 
values  of  both  the  complex  wave  number  and  characteristic  impedance  determined 
from  the  new  method  compared  well  with  classical  methods;  the  values  obtained  by 
the  various  methods  agree  with  each  other  to  within  26%. 


Figure  4.4  |Zcharl/Poc  vs  ^  For  AU-  Methods 


Figure  4.5  zchar  vs  f  For  All  Methods 


APPENDIX  A 

.'1ETH0D  OF  SOLUTION  OF  THE  DISPERSION  RELATION  AND  ITS  IMPLEMENTATION 

The  mehtod  used  to  find  the  roots  of  the  dispersion  relation  (Eq.  (2.51) ) 
was  Ward's  method  [Ref.  19,  pp.  169-170].  This  method  is  not  the  only  method 
available  to  find  the  roots  of  a  complex  transcendental  equation  but  is  ^asy  to 
understand  and  to  implement,  as  no  derivatives  are  needed. 

1.  WARD’S  METHOD 

If  we  let  the  polynomial  f(z)  be  a  function  of  the  complex  variable 
z  =  x  +  iy  such  that 

f(z)  =  u(x,y)  +  iv(x,y)  ,  (A. 1.1) 

it  can  be  shown  that  f(z)  has  a  zero  a  *  b  +  id  if  and  only  if  u(b,d)  =  v(b,d)  = 
0.  This  pair  of  simultaneous  equations  in  u(x,y)  and  v(x,y)  may  be 
conveniently  solved  in  two  ways.  First,  one  can  seek,  to  minimize 

u2(x,y)  +  v2(x,y)  (A.1.2) 

or,  secondly,  on  can  minimize 

w(x,y)  =  |u(x,y)|  +  |v(x,y)|.  (A. 1.3) 

Since  no  derivitives  are  employed,  this  method  seems  well  suited  for  finding 
multiple  roots. 

Ward’s  technique  is  based  on  the  knowledge  that  in  every  neighborhood  of  a 
point  (x,y)  such  that  w(x,y)  *  0,  there  exists  a  point  w(x’,y’)  such  that  w(x', 
y')  <  w(x,y).  At  a  given  step  in  the  process,  which  is  iterative,  the  current 
value  of  w(x,y)  is  oompared  with  the  vlues  w(x+h,y),  w(x-h,y),  w(x,y+h)  and 
w(x,y-h)  until  a  smaller  value  is  found  (if  possible).  The  coordinates  yielding 
this  snaller  value  of  x  and  y  are  used  in  the  next  iterative  step.  Should  these 
four  points  yield  no  decrease  in  w  =  |u|  +  |v|,  the  increment  h  is  decreased 
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by  some  amount  and  the  process  is  repeated.  As  with  other  methods  which  can  show 
convergence  within  certain  limits,  there  is  no  guarantee  of  convergence 
with  this  method. 

2.  IMPLEMENTING  WARD’S  METHOD 

In  order  to  apply  Ward's  method  to  find  the  roots  of  the  dispersion 
relation,  it  must  be  separated  into  its  real  and  imaginary  parts.  To  do  so  note 
that  the  parameters  kjz  and  k2z,  and  at(ui)  are  functions  of  w  ,  they  may  be 
separated  into  their  respective  real  and  imaginary  components  (these  three  terns 
also  have  complex  parameters  buried  within  them  which  are  freguency  dependent). 
Although  tedious,  it  is  necessary  with  this  method. 

Once  each  parameter  has  been  reduced  to  its  simplest  complex  form  (real 
and  imaginary  portions),  these  can  be  inserted  back  into  the  dispersion  equation 
and  by  use  of  complex  arithmetic,  the  dispersion  relation  may  be  recast  in  the 
form 

Re  {Dispersion  Relation}  +  jIm{Disperson  Relation}  =  0 
and  Ward's  Method  of  Eq.  (A. 1.3)  may  be  new  implemented  in  a  computer  algorithm. 

The  numerical  analysis  program  written  using  Ward's  Method  required  the 
following  inputs  for  each  sample  thickness: 

(1)  Measured  Quantities 

a)  fc  (measured  in  the  cavity) 

b)  0  (measured  in  the  cavity) 

(2)  Cavity  Parameters 

a)  Length  of  the  cavity,  L.  (in  x-direction) 

b)  Width  of  porous  material,  d2»  (sample  thickness) 

c)  Width  of  remaining  bulk  fluid,  dj. 

d)  Mode  number  in  the  long  direction,  (n  ■  0  for  the  (010]  mode, 
n  *  1  for  the  (100)  mode} 
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Fluid  and  Porous  Solid  Parameters 


(3) 

a)  Adiabatic  fluid  sound  speed,  Cf.  (®  23*C,  345.2  m/s) 

b)  Fluid  bulk  density,  pf.  <@  20*C,  1.21  Kq/m3) 

c)  Adiabatic  ratio  of  specific  heat  of  fluid,  y.  (<?  20*C,  1.402) 

d)  Prandtl  number,  Pr.  20*C,  0.733) 

e)  Porosity,  P.  ( -  1.0) 

(4)  Program  Parameters 

a)  Angular  frequency  increment,  Aw. 

b)  Dampinq  frequency  increment,  AS. 

c)  An  endinq  criterion  for  w(w,  $). 

(5)  Initial  Estimate  of  Rn^ 

a)  If  no  a  priori  knowledge  of  material's  flow  resistance,  then 
assume  very  hiogh  or  low  value  as  the  starting  point. 

b)  Perform  a  DC  flow  resistance  experiment  to  extract  a  mean  value 
for  Rpc* 

Once  these  quantities  are  entered,  the  proqram  calculates  its  first  w(u»,8)  based 
on  the  initial  inputs  of  fc  and  0.  This  value  of  w(w, 8)  will  be  compared  to  the 
four  incremental  values  of  w(w  ±  Aw,  S  t  AS),  that  is,  it  compares  the 
initial  w(  w, 6)  with 

w(  w*- Aw,  S) 
w(  w-Aw,  6) 
w(  w,  S+AS) 
w(  W,  0r- AS) 

Which  ever  is  the  lowest,  those  coordinates  in  w,S  space  are  chosen  as  the  new 
values  and  the  process  is  continued  until  w(w, 6)  is  less  than  the  endinq 
criterion.  The  output  converts  the  roots  of  w  and  6  to  fc  and  0. 
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This  process  is  continued  for  each  sample  thickness,  and  the  resulting 
values  can  be  plotted  either  as  Q  vs  fc  or  as  fc  vs  sample  thickness  and 
compared  with  the  respective  curve  after  loss  corrections  to  Q  are  made.  An 
evaluation  is  made  as  to  whether  Rqc  is  too  hiqh  or  too  low,  the  value  of  Rqq 
changed  and  the  process  repeated  for  this  new  value  of  R^c  f°r  all  samples. 
This  iterative  process  is  completed  when  the  roots  are  found  for  a  value  of  Rq£ 
which  provides  the  best  fit  curve  through  the  measured  data.  With  this  value  of 
Rqq  and  the  roots  for  each  sample  thickness  determined,  the  characteristic 
impedance  and  complex  wave  number  may  be  calculated  by  use  of  Eg.  (2.56)  and  Eq. 
(2.57). 


123 


APPENDIX  B 


DETERMINATION  OF  TOE  THEORETICAL  MODAL  PRESSURE  DISTRIBUTION  IN  A  PLANE 
VJAVE  RESONATOR  WITH  AN  ACOUSTICALLY  ABSORBING  MATERIAL 


Jut 


The  equation  that  governs  the  pressure  distributuion  in  the  bulk  fluid 
(region  1)  is:  j— 

g~  =  -jolpfcos[kxxl  £10  cos(klz(z-d1)] 

The  equation  that  governs  the  pressure  distribution  in  the  porous  material 
(region  2)  is,  after  substituting  Eq.  (2.62)  for  #20: 


(B.1.1) 


*P2 


cos[k  djJ  ^ 

-j“PfCoslkxx]^i0  cSsT^TT  costk2z(z+d2)] 

2z  2 


Jut 


(B.1.2) 


With  the  value  of  R^  known,  along  with  the  roots  of  the  dispersion  relation  for 
a  given  sample  thickness,  the  pressure  distribution  in  region  1  and  region  2 
can  be  calculated. 

In  the  experimental  measurement  of  the  modal  pressure  field,  the  quadrature 
component  of  the  pressure  was  zeroed  at  the  reference  position  of  (0,0,1)  to  set 
a  reference  point  such  that  both  theory  and  experiment  could  be  compared.  The 
0,0  indicates  the  grid  point  (0,0)  and  the  1  indicates  that  the  receiver  probe 
tube  was  located  in  the  ceiling  of  the  cavity.  To  set  the  reference  at  x  =  0,  z 
=  dj  such  that  the  phase  of  the  pressure  distribution  at  this  point  is  zero, 
requires  that 

$10  =  1  /-ATO(-ay/S)  .  (B.1.3) 

To  calculate  the  theoretical  modal  pressure  distribution,  start  at  the 
reference  location  of  x  =0  and  z  =  d^  and  incrementally  move  away  from  this 
position  in  increments  of  Ax  and  Az  such  that  the  theoretical  grid  points  match 
those  existing  in  the  actual  cavity. 
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